PROBLÈMES CORRIGÉS 


Il n’est pas de problème qu'une absence de 
solution ne finisse par résoudre. 


Aphorisme attribué à Henri QUEUILLE 


PROBLÈME N° 1 : Coefficient d’échange 
Énoncé 


Le transfert de chaleur entre deux fluides s'effectue à travers un tube d'acier de 
diamètres intérieur/extérieur 18 / 21 mm. 
On donne : 


- côté intérieur : hy = 1000 W / m°K ; température moyenne de mélange T; = 10°C 
- côté extérieur : hy = 2000 W/m°K ; température T, = 25°C 
- acier: À=46W/m.K 


1. Calculer le coefficient global d'échange k. 
2. Après un an de fonctionnement, on estime avoir une résistance d'encrassement 


R,=4.10 1 W~! m? K. Déterminer le nouveau coefficient d'échange global. 


3. En attribuant une efficacité de 1 au tube neuf, que devient cette efficacité au bout 
dunan? 
4. Quel est alors le flux échangé dans un tube de longueur L = 1m? 


Solution 


1. 
La paroi du tube a pour épaisseur : 


SË 1,5 mm — 15.107 m 


Comme elle est mince par rapport aux diamètres, en négligeant sa courbure on peut 
calculer k à partir de la formule (6.2a) relative à une paroi plane, avec R, = 0 : 
DD ee ed tr MS il 

+—+—= + + 
k hy À h, 1000 46 2000 


e 








nn 0,03.10 *+0,5.10 *=1,53.10 7° 


k 
On constate que la résistance thermique e/ À de la paroi ne représente ici que 2% de 


la résistance totale. Enfin : 





k=653W/m° K 











2. 
En présence d’une résistance d’encrassement, on applique maintenant la formule (6.2a) 
complète : 
1 NI 1 
—=— + R+ nn 
k hi A hy 


=(1+0,4+0,03+0,5)107°=1,93.10 








k=5I8W/m°K 








3. 
L'efficacité dont il est question dans l’énoncé doit être comprise comme un rapport 


D cel / Pmax (définitions 3.5 et 4.37), soit ici : 
T Pran = kian 518 

Pneu K 653 
E = 0,793 


E 





neuf 














4. 
La surface d’échange Z n’est pas la même des deux côtés. Suite à la question 1, on 
calcule une valeur approchée de 4 par (6.10) ($ 6.2.3 9): 


Lex E 3, ] 








TIT 


PSGTIO m? 
Le flux échangé est donné par (6.2c) qui s’écrit avec les notations de l’énoncé : 


D=kZ(T-T;,)=518X61.10 °x15 





PD -475V 











Commentaires 


Cet exercice très élémentaire a surtout pour objet de matérialiser des ordres de 


grandeur. 
Pour les étudiants, sa principale difficulté réside dans le calcul de l’épaisseur de la 


paroi, qui n’est pas D — d comme on le lit souvent !! 


PROBLÈME N° 2 : Méthode NUT 
Énoncé 


Un échangeur à contre-courant fonctionne dans les conditions suivantes : 


T.e = 350°C T po = 120°C 
T., = 200°C Ts = 290°C 
Qimin = dif Puissance P = 415 kW 


1. Quelle est la puissance ëchangëe si on fait travailler l’échangeur en mode co- 
courant, avec les mêmes températures d’entrée et les mêmes débits ? (Utiliser la méthode 
NUT). 

2. Quelles sont les nouvelles tempëratures de sortie 2 


Solution 


1. 
En contre-courant, avec Gjin = qıf , On a pour efficacité (cf. 3.11) : 


Ti” Tje 290 — 120 


E = — — 0,74 
Ty — Tje 350 — 120 
R= Tee Tos Z 350 — 200 = 0,882 
T-T 290 — 120 
et, à partir du tableau 3.1 : 
NUT = — Pan e 





1- R l-E 
En circulation co-courant, les débits n'ëtant pas modifiés, les coefficients d’échange ne 
le sont pas non plus. On garde donc le même NUT (vu que NUT =k 7 qy, ). Par contre, 


la nouvelle efficacité E” s’écrit (tableau 3.1) : 


pe U-(1-R) NUT|} 


I+R 
Il vient, après remplacement de R et NUT par leur valeur : 
E' = 0,526 


Puisque les conditions d'entrëe sont identiques dans les deux cas, la nouvelle 
puissance ®' est telle que : 








Per 
D E 
ceci d’après la relation (3.12). Alors : 
D'=415 2e 
0,74 
D' = 295 kW 











c’est-à-dire 70% de la puissance en contre-courant. 


2. 


La nouvelle température de sortie froide Th s’obtient à partir de la nouvelle efficacité : 


' 


NË Th më T e 
Në Tje 


' 


Th =T p+ E (To The) 
= 120 + 0,526 (350 — 120) 





G 


TeC 











et la nouvelle température de sortie chaude Tos 


T -T 














R = K—S = 0,882 inchangé 
Tp- Th 
Tis =T- R(T- Tje) 
= 350 — 0,882 (241 — 120) 
TT 243,3°C 
Commentaires 


à partir de R : 


II y a d’autres façons de résoudre ce petit exercice, par exemple en utilisant les 
résultats du chapitre 2. Mais la démarche préconisée permet de se roder à la méthode NUT 


dans un cas simple. 


Avec cette valeur de R, on est déjà dans la zone asymptotique d’un échangeur co- 
courant (fig. 3.1), et les températures de sortie sont donc très voisines. 


PROBLÈME N° 3 : Échangeur bitube 
Énoncé 


Pour refroidir un débit de 9,4 kg/h d'air de 616 “Ca 178 °C, on le fait passer dans le 
tube central d’un échangeur bitube à contre-courant de 1,5 m de long, de 2 cm de diamètre et 
de faible épaisseur. 

1. Calculer la puissance calorifique à évacuer. On donne pour lair: 
Cpe= 1060 J/kg K. 


2. Le fluide de refroidissement est de l’eau, qui pénètre dans la section annulaire à la 
température de 16 °C avec un débit de 0,6 l/mn. Calculer la température de cette eau à la 
sortie de l'échangeur. On prendra Cpp = 4180 J/kg K. 

3. Calculer le coefficient d'échange h, côté fluide chaud (on ne tiendra pas compte 
d'une éventuelle correction en 4/ Hp). 

4. Déterminer l'efficacité de cet échangeur, puis son NUT. En déduire le coefficient 
d'échange global, puis le coefficient d'échange hy côté fluide froid. 

5. La paroi extérieure de l'échangeur est isolée. Quelle est approximativement 
l'épaisseur b de l’espace annulaire qui permettrait d'obtenir cette valeur de hy 7 (On 


admettra d’abord l'écoulement laminaire, et on vérifiera ensuite cette propriété). 
Solution 


1. 
Le flux total peut se calculer côté chaud : 
DSgAI ST 
D’après les données, le débit thermique unitaire chaud est : 
den Cpe = t x 1060 
4 =2,77 W/K 
et alors : 
D = 2,77 (616 — 178) 








p = 1213W 











2. 
Le calcul du flux total côté froid va maintenant nous donner T p . 


P=qy(Ts Tf) avecici: 





0,6 
If = Amf Cpf = x 4180 (amf = 0,3 kg/mn) 
dy =41,8 W/K 
d’où 
D 1213 
Te ZTy, + — = 16 + — = 16 + 29 
PU 41,8 


dif , 





Ti=4 °C 











3. 
Il faut d’abord connaître le régime d’écoulement de l’air, donc le Reynolds côté chaud. 
La température moyenne de l’air est approximativement ($ 6.2.1) : 
Tet Ts _ 616 + 178 
2 
A cette température, les tables donnent : 
pe =0,525 kg/m? ; v, =6,20.10 °m?/s 


d RX(2A10 AY 


<T, >= = 397°C = 670 K 








La section du tube est : $, = 7 -— = —————— — 3,14.10 i m? 
On en déduit la vitesse débitante : 
_ am _ 94 1 
E Pë 5. 3600 0,525 431410“ 
V.=15,8 m/s 
d’où le nombre de Reynolds : 
V.d 15,8x2.10 7? 
Pe eee = 5100 


c 


Ve 6,20.10”5 
II s'agit d'un rëgime de transition. On peut donc utiliser la formule (4.26), en notant 
que le rapport d7 L = 27150 est négligeable (I'ënoncë nous demande également d'ignorer la 


correction en UJ 4p): 
= CHÓ (Re 275 — 125) pr 213 
Re C C 


Pour calculer h,, il est un peu plus rapide ici de passer par le nombre de Nusselt 
(4.10d, $ 4.1.5) : 


Nu, St, Re, Pr,— 0,116 (Re2/ 3 125) Pr!” 
A 670 K, le nombre de Prandtl de l’air est : Pr, = 0,68 . On trouve : 


St, 


c 





Nu, 17,25 

Toujours à 670 K, la conductivité de l’air est : A, = 0,0505 VImK. 
h.d 17,25 x 0,0505 

AE NE unes 
À 0,02 





h.=43,5W/m° K 











4. 

On constate que : 

Qt min = dre = 2,77 W / K , d’où l'efficacité (formule 3.10) : 

E = Te Lu 616 — 178 438 
1212.  6l6-16 600 








E = 0,73 











D’après le tableau 3.1, pour un échangeur à contre-courant : 
nes ere 
1-R LE 
et dans le cas présent : 
VA. dt min pë 2,77 








R= = 0,066 
dr max 41,8 
donc : 
NUT = 1n 17 9:066 x 0,73 = 1,07 Ln3,52 
1 — 0,066 1— 0,73 
NUT = 1,35 











De la définition du NUT (3.14a) on tire alors : 
NUT qy min 


pa 
Puisque l’épaisseur du tube central est faible, on ne fait pas la distinction entre surface 
d’échange côté chaud et côté froid, et on néglige la résistance thermique de la paroi. Donc : 
Z = surface latérale du tube = xd L 
Z=xx2.10 ?x1,5 = 0,094 m 
_ 1,35 x 2,77 
0,094 


k=398W/m° K 


pe 














Le coefficient d'ëchange global s'exprime aussi à partir de (6.2a) (e et R, étant 


nëgligës) : 
Dr 1 1 1 1 





+ — = + —= 
k h h f 43,5 h f 39,8 
On en déduit : 





hp = 500 W/m° K 











5. 
Dans un échange eau-air, la température de paroi est proche de celle de l’eau. Dans le 


cas présent, celle-ci varie peu. On peut donc admettre la condition T,= cte. La paroi 
extérieure (concave) étant isolée, et lécoulement supposé laminaire, la formule (4.44) 
s'applique. Elle donne Nu en fonction de R, / R}. 

D'autre part, Nu = hp Dp / Åp, avec ici D,=2(R,-—R;)=2b (formule 4.42). En 
partant d’une valeur arbitraire mais raisonnable de b, on peut par approximations successives 
ajuster Nu avec la valeur donnée par (4.44). 


La température moyenne approchée de l’eau (fluide froid) est : 


Tat IT Kë 16 + 45 


<T; > = = 30,5°C 


La conductivité correspondante (voir tables) est À; = 0,612W/mK. 


aus. Ro _10+3 
Essayons avec R, — R; =b = 3 mm (ce qui fait —2 = 
R; 10 





=1,3). 


500x2x 3.107 
0,612 
D’après (4.44), pour obtenir cette valeur de Nu, il faudrait un rapport R» / R; = 1,05 


environ, soit b = 0,5 mm, ce qui est trop faible. Essayons avec b un peu plus élevé, pour 
augmenter Nu. 


Par exemple, avec b = 3,2 mm (soit R, / R; = 1,32) , Nu = 5,23 ce qui correspond 
à peu près à un rapport R/R; =1,35 dans (4.44). On admettra donc comme valeur 
approchée : 


Nu = 4,92 








b = 3,2 mm 











On vérifie enfin le Reynolds : 
Vf Dr _ amf 2b 
ye SË 


avec V=0,083.107 RE D T;>=30,5°C et, en première approximation 





Reç — 


Sf =bxz(d+2(b/2)) surface d’un rectangle de hauteur b et de longueur égale à la 


circonférence moyenne de l’annulaire, soit z( d + b) (N.B. le diamètre moyen est d + b, et 
le diamètre extérieur d + 2b, voir Problème 1). 





0,6 
D'autre part, g,yç = 0,6kg / mn = F kg/s 





amf 2b 1 0,6 2 1 
Rep = = 3 —3 55 
ps br(d+b)v; 60ox10 (20 t3,2)10 0,083.10 
Re ; = 331 


L’écoulement est bien laminaire. 
Commentaires 
Dans cet exercice, on doit en particulier chercher une caractéristique géométrique de 


l’échangeur permettant de respecter les conditions thermiques imposées. 
On a en plus l’occasion d’aborder régime de transition et écoulement annulaire. 


PROBLÈME N° 4 : Cheminée 
Énoncé 


On veut estimer la chute de température des fumées dans une cheminée, en 
considérant le conduit comme un échangeur dont les fumées constituent le fluide chaud, et 
lair ambiant le fluide froid. On admet que la température T, de l’air est constante le long de 


la paroi extérieure de la cheminée. On désigne par k le coefficient global d'échange à travers 
la paroi. 


1. En adaptant le calcul d’un échangeur co-courant au cas particulier ci-dessus 
(T; =T,= cte ), montrer que la température des fumées dans la cheminée obéit à la loi : 


T.—T 
që sen Es) 
T T dte 


ce ta 
2. Le conduit est cylindrique, de diamètre D et de longueur L. Écrire la température 
T,, de sortie des fumées. 


3. Calculer T,, avec les valeurs suivantes : L = 20m ; D=30cm; T,= 320°C ; 


T, = 10°C ; k = 20 W /m°K. Pour les fumées © Ame = 0,5 kg/s ; Che =1050 J/kgK. 


Solution 


1. 
A travers un élément de paroi dS de la cheminée s'ëchappe un flux de chaleur : 
dp =k(T.-T; )dS =k (T, -T ) dS 


d est aussi la chaleur perdue par le fluide chaud : 
dB =- dte dT, 


On regroupe les deux équations : 
dT, 2E k dS 
TT, dte 
Nous avons ici T, = cte , de sorte que : 
dT,=d(T,—T,) 
d’où en intégrant : 


T.—T,= cte x e|- La s) 





dte 
A entrée, S =0 et T, —T,,, donc cte = T e — T}, €t 
T—T k 
——C = exp - — s) C.Q.F.D. 
Tee = Ta dte 
2. 


A la sortie, T, = T,, et S = Z (surface latérale totale de la cheminée) : 


ZË -xDL 
La formule de la question 1 devient : 


k 
Tes = Ti + Tee- T ) e|- + 27 
dtc 


3. 
Les calculs numériques donnent : 
Gje = Ame C pe = 0,5 X 1050 = 525 W / K 


Ta = 10 + (320 — 10) exp- À z X0.3 x 20) 








T.. = 161 °C 








La chute de température est donc de 760 °C environ dans la cheminée. 


Commentaires 


1. La formule donnant T,, (question 2) est recommandée par le Centre Scientifique et 
Technique du Bâtiment (CSTB) pour le calcul des cheminées. 

2. Dans la question 1, on ne peut pas écrire le bilan local sur le fluide froid 
(dB = Gif dT; ) du fait que dT; = dT ,=0 (T,= cte ), ce qui est équivalent à qfy >. On 
se retrouve dans la même situation qu’avec un évaporateur ou un condenseur ($ 2.4 et 3.5). 

3. L'hypothèse T,= cte est acceptable (T, n’est pas ici une température de mélange 
mais la température du fluide extérieur loin de la paroi) ; elle conduit à une formule simple 
qui donne des estimations numériques correctes. 

4. Dans une revue technique, la formule de la question 2 était accompagnée du 
commentaire suivant : X est la section de la cheminée !! Ce qui donnerait ici une chute de 
température inférieure à 1 °C. De l’importance de la rigueur dans le vocabulaire !! 


PROBLÈME N°5: Échangeur à changement de phase 
Énoncé 


De l’eau froide circule dans un tube de chaudière à condensation. Sa température 
d'entrée est Ty SI8“C et son débit q,yç = 400 kg/h. Le réchauffage est assuré par 
condensation de vapeur d’eau à l'extérieur du tube, à la température T,= 104 °C. 

On donne: diamètre intérieur d = 12,5 mm ; diamètre extérieur D = 16 mm , 


longueur L = 2,4 m ; conductivité de la paroi Ap =46 W/mK. Pour l’eau, on admettra 


dans la gamme de température considérée: Vy =0, 710 mis ; Pr; =5,5 ; 
Cp = 4180 J/kgK. 


1. Calculer le coefficient d'échange hy à l'intérieur du tube. 


2. On donne le coefficient h côté vapeur : h, = 8000 W /m° K. Calculer le coefficient 
global d'échange k. 

3. Calculer le NUT et l'efficacité de l’appareil. 

4. Déterminer la température de sortie d’eau Tp, puis la quantité de chaleur Q 


récupérée annuellement grâce au dispositif à condensation, si l’on considère que la saison de 
chauffe dure 150 jours et que la chaudière fonctionne 5 heures par jour. 


Solution 


NË 
II faut dëja savoir quel est le rëgime d'ëcoulement. Pour cela, calculons Re : 
e, 





Re 
V Sv zxzdv 
z3 
Ici, T A eu RS m? /s 
3600 
0,4 
Re = 4 x —— x e z 
3600 mx12,5x10™x0,7 x107 
Re = 16170 


Avec cette valeur de Re, le régime est turbulent et hydrauliquement lisse. Donc, St est 
donné par (4.21). On a bien Pr > 0,66 . D'autre part, L/d = 2400/12,5 >60 . Les 
conditions de validité de la formule sont donc satisfaites. 


Pour le fluide froid : 


St ; = 0,023 Re” 92 Pr” 95 = 0,023 x16170 7°? x 5,57 06 


h 
Sip LI910”5 et aussi Set d'où : 
pPCpV 





4 q 
h,=St,pC,V=St;pcC Z 
Í f P f P xd? 





Ona p= 10? kg/m? pour l’eau et d = 12,5 mm, ce qui donne : 











hş = 4500 W/m°K 





2. 

L'ëpaisseur de la paroi est : 

= 2 16 —12,5 — 175 mm 

La question ne précise pas de quel côté on demande k. Cela sous-entend que l’on 


raisonne comme avec une paroi plane ($ 6.2.1. et 6.2.3 ® ), et k s'ëcrit : 
I1 1l e l 1 1751075 1 
= —+—+ + + 


= = = 3,85.10 71 
k h À hy 4500 46 8000 





C 





k = 3000 W/m°K 











3. 
L’échangeur est un appareil à fluide isotherme, pour lequel R = 0 ($ 3.5.2), avec ici 


Qi min = lef - Ainsi : 





NUT = KZ/qy 
où Z est la surface moyenne d’échange ($ 6.2.2) : 
pepe L= gx EEE x 10 7x 2,4 = 0,107 m? 


Le débit thermique unitaire de l’eau a pour valeur : 


400 


etle NUT vaut : 





NUT = 0,69 











L'efficacité se trouve sur le tableau 3.1. Pour R=0 , quelque soit le modèle 


d’échangeur : 
E = 1 — exp(— NUT ) = 1 — exp(— 0,69) 














Toujours pour la raison que q; min = qy» E est l'efficacité côté froid (formule 3.11), 
avec Te =T, = cte : 
pb 
ST 


On en tire Ty, : 
Ty = 61°C 














La puissance thermique du système s'exprime par le bilan sur le fluide froid : 
P = qy (Ts - Ty )= 465 (61 — 18 ) = 20000 W 
D = 20 kW 


Sur une saison de chauffe qui compte 150 x 5 heures, on peut donc récupérer par 
condensation au maximum : 
Q = 20 x 150 x 5 kWh 








Q = 15000 kWh 








Commentaires 


Il s’agit encore d’un problème allégé où on ne fait mention ni des températures de 
référence, ni de la distinction entre coefficient d’échange global côté froid et côté chaud, et où 
le calcul de h côté vapeur n’est pas demandé. 

La température de condensation de 704 °C sous-entend que la vapeur est en légère 
surpression par rapport à la pression atmosphérique. 

Les questions 3 et 4 peuvent aussi être traitées avec les résultats relatifs aux 
échangeurs à fluide isotherme (§ 2.4). 


PROBLÈME N° 6 : Échangeur à faisceau de tubes et calandre 
Énoncé 


Dans la sous-station de chauffage collectif d’un immeuble on désire installer un 
échangeur à faisceau tubulaire et calandre, destiné à porter de 40 à 60 °C un débit d’eau de 
20000 kg/h. Le fluide primaire qui circule dans les tubes est de l’eau surchauffée arrivant à 
180 °C, à raison de 10000 kg/h. Les tubes ont un diamètre intérieur d = 20 mm ; la vitesse 
d'écoulement adoptée est telle que Re = 10000. Le coefficient d'échange global k est estimé à 


450 W/m°K. On admet pour l’eau surchauffée les caractéristiques thermophysiques 
suivantes : C,.= 4315 J/kg K ; p = 920 kg/m? ; H= 19.10 kg/m.s. 


1. Calculer la puissance P ëchangëe et la température de sortie du fluide chaud. 
2. L'ëchangeur est à contre-courant, avec une seule passe sur chaque fluide, les tubes 
étant montés en parallèle (modèle de la fig. 1.5 sans les chicanes). Déterminer : 
- la surface d'échange X nécessaire 
- la vitesse dans les tubes 
- la section totale des tubes 
- le nombre de tubes et la longueur du faisceau 


Solution 


1. 
Calculons tout d’abord les débits thermiques unitaires. 
- Sur le fluide froid : 


If = amf Cpf 
On prend généralement pour l’eau Cpe = 4180 J/kg K (voir tables en fin de volume). 


20000 
z AN sh 
If = -3600 


qış = 23200 W / K 


- Sur le fluide chaud : 
10000 
Gje = Ame Cpe = 3600 


qı = 12000 W / K 





x 4315 


La puissance thermique s’obtient à partir des données relatives au fluide froid : 
D = qy (Ty — T fe ) = 23200 X(60 — 40 ) = 46400 W 
D = 464 kW 














On tire maintenant 7,. du bilan sur le fluide chaud : 
B = qie (Tee ~ Tes ) 

464000 = 12000 x ( 180 —T,.) 

T., = 180 — 38,7 

T. = 141,3 °C 














2. 


se Les calculs précédents montrent que nous avons ici : 





dt min — itc 
et l'efficacité s’écrit : 
_ Te Tes _ 38,7 
Te Te 140 
E = 0,276 
Quant au facteur de déséquilibre : 
Gr. _ 12000 
qy 23200 
R=0,517 


Du tableau 3.1 on tire le NUT : 





NUT = 1 Ln I-RE 10 
1-R 1-E 1-0,517 
NUT = 0,35 
et aussi : 
NUT =k£Z/q, 


ce qui nous donne la surface d’échange : 
s- 0,35 x 12000 


450 
S =9,33 m 














© Le nombre de Reynolds est imposé : 
pVd 
u 


Re = 10000 = 


d’où la vitesse V : 
pa 19.10 7° x 10* 
920 x 20.107? 





V =0,103 m/s 











1 — 0,517 x 0,276 


1 — 0,276 


V La vitesse V est la même dans tous les tubes. Si q, désigne le dëbit-volume total, 


la section totale des tubes est donc : 
dy _ Ame 1 _ 10000 1 1 
V p V 3600 920 0,103 


S = 2,9.107? m? 

















À Le nombre de tubes nécessaires pour assurer le débit demandé est égal à : 


—Z 
ga OO A 


zd? /4 xx(2.10 ? } 
Comme il faut bien avoir un nombre entier de tubes, on prendra l’entier 
immédiatement supérieur : 








N = 93 tubes 











Si la longueur du faisceau est L (c’est donc aussi la longueur de chaque tube), la 
surface totale d'échange a pour valeur : 








EzNxdL 
On en dëduit : 

E — 933 

= Nad 93x7 x 0,02 
L=1,60m 











soit une longueur totale de tubes de 749 m environ. 
Commentaires 


1. L'exercice aborde d’une manière encore élémentaire le dimensionnement d’un 
échangeur à faisceau de tubes. 

2. Les étudiants qui n’ont jamais vu un échangeur tubulaire n’imaginent pas toujours 
très bien la structure de l’ensemble et les propriétés qui découlent du montage des tubes en 
parallèle (FEMM, $ 7.6). 

3. Le NUT de l’échangeur est faible. Le modèle de l’appareil n'a donc pas beaucoup 
d'importance ($ 3.4.3 et fig. 3.1). 

4. Un calcul approché de X pourrait être effectué à partir du bilan global en 
linëarisant la distribution de température de l’eau ($ 6.2.1) : 

TetTs 
2 

C’est plus rapide et on trouve ici pratiquement le même résultat (parce que NUT est 
petit). Mais ce n’est pas toujours comme ça. Alors prudence !! 

5. Les paramètres k et Z se réfèrent implicitement à la surface intérieure des tubes 
puisqu'on ne donne pas le diamètre extérieur dans l’ énoncé. 


B=k£(T,-<T; >) avec <T; >= 


PROBLÈME N°7: Échangeur à plaques 
Énoncé 


De l’air chaud à 80°C produit par une installation industrielle est évacué à l’extérieur 
à raison de 5000 kg/h. Il traverse préalablement un échangeur servant à préchauffer à 20°C 
lair neuf admis dans le bâtiment à une température d'entrëe de 5°C. 

L'ëchangeur est un appareil à plaques planes et à courants croisés, sans brassage des 


fluides, dont la surface d'échange est X = 50 m”. Dans la plage de fonctionnement prévue, 


son coefficient global d'échange k est estimé par le constructeur à : 
ks 3,6 qm + 65 


avec qy en kg/s etken W/m°K. 


1. Estimer a priori sur quel circuit le débit thermique unitaire est minimal. 

2 En procédant par approximations successives, déterminer quel débit d'air neuf cet 
échangeur peut porter à 20°C. 

3. Calculer la température de sortie du fluide chaud et la puissance de l'échangeur. 

P.J. Abaque E = A(NUT), fig. P.7. 


0,6 





E A 
1 2 3 4 5 NU 


T 


FIG. P.7. — Echangeur à courants croisés, fluides non brassés 


Solution 


1. 

Ona AT =(20 — 5) = 15 °C sur le fluide froid. Le fluide chaud entre à T.e = 80°C , 
ce qui doit permettre un AT nettement plus important. Les fluides étant les mêmes des deux 
côtés, qms est certainement supérieur à qme» d'où Gy > qre Car la chaleur massique C, 


dépend très peu de la température. Par conséquent, a priori : 


Qi min — Are 
puisque le fluide qui accuse le plus grand AT est celui qui a le plus petit débit thermique 
unitaire ($ 3.1). 


2. 
Les données ne portent que sur trois températures aux bornes et un débit. Le débit 
inconnu gjy ne peut pas être calculé directement ; l'efficacité non plus puisque ici E = E, et 


on ne connaît pas 7... On va donc procéder par approximations successives à partir d’une 
estimation empirique de g,, en calculant NUT, puis E au moyen de l’abaque ( 1! dernier 
rappel : abaque est un nom masculin, on dit un abaque, fermez la parenthèse) et enfin Ty, que 


l’on va comparer à la valeur imposée 20 °C. 


se Choisissons comme valeur de départ : 
amf = 15000 kg /h = 4,166 kg/s 
None, jé, Sen 
dt min dte 
me = 3000 / 3600 = 1,39 kg / s 
D’après les tables numériques, C „e = 1006 J/kg K très peu dépendant de T. Soit : 


Gje = Ame C pe = 1400 W / K 





Dans ces conditions : 
k =3,6 qmf + 65 = 3,6 X 4,166 + 65 = SON Im” K 
80x50 


1400 
R = ic = mc 





, 


puisque C e= Cr 
qa An pe pf 
R = 1,39/ 4,166 = 0,333 


En reportant sur l’ abaque, on trouve pour l'efficacité: E = 0,85 


Puisque q; min = Gj, NOUS avons (formules 3.10) : 


d’où l’on tire : 
Th =Tp+ RE (Tee Tre) 
= 5 +0,333 X 0,85 x(80 —-5) 
soit : 
Ty, = 26,25 °C 
Cette valeur est supérieure à la température imposée (20 °C) ; on a donc pris un débit 
qmr trop faible. 


9 Recommençons l’opération avec : 

Gj = 25000kg / h = 6,945 kg / s 

La trame du raisonnement est identique ; q,, est inchangée = 1400 W / K 
90 x 50 _ 32 


k = 3,6 x 6,945 +65=90W/m° K : NUT = 
1400 





R=1,39/6,945 = 0,2 
et l’abaque donne : E = 0,9 


T=5+0,9Xx0,2X (805) 
Ty = 18,5 °C 
Nous nous sommes nettement rapprochés de 20 °C, mais maintenant Ty, est un peu 


trop faible ; donc q,r est un peu trop grand. 


V Faisons une troisième tentative avec Inf = 22500 kg/h =6,25 kg/s 


k — 3,6 X6,25 +65 -87,5 W/m? K 
NUT = 87,5 x 50 = 3,125 
1400 


R=1,39/6,25 = 0,222 
Sur l’ abaque, on observe que E est pratiquement inchangée = 0,9 
T=5+0,9x0,222x(80 —-5) 


Ty, = 20 °C : La convergence est atteinte. La valeur demandée est donc : 





Amf = 22500 kg/h = 6,25 kg / s 











3. 
Nous pouvons par exemple déduire 7,, du bilan global de l’échangeur (3.1) : 
dte (Tee TE Tes )= dif (Th E T jo ) d’où : 
6,25 
Tes EE E i 
dte 1,39 
puisque qf / Gje = Qmf / me - Finalement : 





(20-5) 








T..= 12,5 °C 








Le bilan sur le fluide chaud (ou froid) donne enfin la puissance Æ de l’échangeur : 
P = qı (Toe — T., ) = 1400 x ( 80 — 12,5) 





PB = 94500 W 











Commentaires 


Cet exercice illustre l’un des rares cas où l’on doive procéder par itération avec la 
méthode NUT . 
Pour le calcul de Ty, (question 2), on peut aussi calculer préalablement T,, à partir de 


E puis déduire Ty, du bilan enthalpique de l’échangeur. 


PROBLÈME N° 8 : ÉCHMEMbNAEMICONEMSOISES 


Énoncé 


Une chaudière à mazout produit des fumées à 360 °C qui sont utilisées pour 
préchauffer de 20 à 120 °C lair nécessaire à la combustion. 

L'installation consomme 1500 kg de mazout par heure. La combustion d’un kilog de 
mazout demande 17,95 kg d'air et produit 18,95 kg de fumées. 

L'appareil utilisé est un échangeur tubulaire à courants croisés, à deux passes côté 
tubes; lair circule dans les tubes, verticalement, et les fumées autour des tubes, 
horizontalement (cette disposition limite les obligations de ramonage). 

La température de paroi T,, des tubes doit être partout supérieure à 160 °C pour 


éviter la corrosion due à la condensation des composés soufrés. 

Les tubes forment un faisceau en ligne, à pas carré, de pas relatif e* = 1,4. Ils ont 
pour diamètre extérieur D = 55 mm et pour diamètre intérieur d = 50 mm. 

Pour chaque écoulement, le nombre de Reynolds est fixé : Re ; = 25000 dans les tubes 


(fluide froid) et Re.= 5000 en calandre (fluide chaud, vitesse de référence = vitesse 


dëbitante en calandre vide ). 
Dans le domaine de température envisagé, on prendra comme caractéristiques 
moyennes des fumées : 


u=3.10 kglm.s  C,=1212J/kgK 
A- 0,033 W/mK p =0,585 kg/m? 


1. Calculer la température de sortie des fumées T,.. 


2. L'ëchangeur peut fonctionner soit avec des entrées opposées, soit avec des entrées 
du même côté (voir figure P8). Quel sens de fonctionnement doit être choisi pour éviter le 
risque de condensation des composés soufrés ? 


3. Calculer à quelles vitesses de référence v? et V? correspondent les nombres de 
Reynolds imposés. 
4. Déterminer le coefficient global d'échange k, côté fluide chaud (négliger la 


résistance thermique de la paroi). 
5. Calculer la puissance ® de l'échangeur, son efficacité E et son NUT. Déterminer 
la surface d'échange X, côté chaud, le nombre de tubes nécessaires et la longueur L de 


chaque tube. 


S e 
Tie Tís 
—agjj—— Re ee a T7 
Të Te Tes ce 
(a) (b) (a) (b) 
FIG. P8 — Échangeur tubulaire à courants croisés. 
Fluide chaud : fumées ; fluide froid : air 
A gauche : entrées opposées. A droite : entrées du même côté 
Solution 
l; 


La température inconnue T,, sera obtenue à partir du bilan enthalpique de 
l échangeur : 

dte (Toe — Tos )= Qi (Tis z T fe ) 

Il faut d’abord calculer les débits thermiques unitaires. 
- Pour le fluide froid (air), le débit-masse est : 

që KE SA 

La chaleur massique de l’air dépend très peu de la température. On l’arrondit 
à Cpr = 1010 J/kg K , soit: 

qyp = Amf Cof = 7,5 X 1010 =7,57.10W/K 
- Pour le fluide chaud (fumées) : 

TRE EE STONGT 

Gje =7,9 x 1212 = 9,57.10° W/K 

On observe au passage que : 

Qt min T Qyg =7,57.10° W/K 


Du bilan enthalpique on tire : 


qy 7,57 
Ts = — (Ts -Ty )= x (120 — 20 ) = 79°C 
cs di ( fs fe) 9 57 ( ) 


tc , 


T, 


ce — 








T,, = 281°C 











2. 

Air et fumées ont des propriétés voisines, et les coefficients d'ëchange à l’extérieur et 
à l’intérieur des tubes doivent être du même ordre de grandeur. Comme température locale de 
paroi, on admettra donc en première approximation : 

T, r (Tair + T fumées )/2 


Avec entrées opposées, on a : 
* côté (a) : air à 20 °C ; fumées à 281 °C, soit : 


20 + 281 
ie E = 150,5 °C < 160°C d’où risque de condensation 
* côté (b) : air à 720 °C ; fumées à 360 °C, soit : 
360 + 120 
He a = 240 °C pas de problème 


Cette solution est à proscrire : il y a un fort risque de condensation vers la sortie des 
fumées. 


Avec entrées du même côté : 
* côté (a) : air à 120 °C ; fumées à 281 °C. Donc : 


120 + 281 
na — = 200°C > 160°C c'est bon 
* côté (b) : air à 20 °C ; fumées à 360 °C : 
20 + 360 
To= T = 190°C > 160C c'est encore bon 


Conclusion : il vaut mieux placer les entrées du même côté 


3. 
Vitesse de référence ve côté fluide froid : 
Il s’agit ici de l’ écoulement de l’air dans les tubes. 





0 
d 
Reç = =25000 donné, avec d = 50 mm 
Vë 
y ç doit être évaluée à la température moyenne de mélange de lair ($ 4.1.5 V), soit : 
Dje 20 + 120 “0 
2 


vr=2,02.107° m°/s 

On en déduit : 

vo 25000 x UE 
50.107 











Vf =]10,1 m/s 








Vitesse de référence ve cotë fluide chaud 
II s’agit maintenant des fumées, sachant que : 


vë D He HO 




















Re, = —— = 5000 avec v, = == 5,13.10 m1s 
ve Pe 0,585 
on en tire : 
y0 Rec Ve. _ 5000 X 5,13.10 7° 
or: 55.10 
VË —4,66 m/s 
4. 


Coefficient d'ëchange h, côté chaud (en calandre) 

Le nombre de Reynolds est donné : Re, = 5000 , avec une vitesse de référence en 
calandre vide, donc compatible avec les formules du parag. 4.3.1. 

Le faisceau est en ligne ; de ce fait il résulte que : ej = ej = e" = 1,4. En admettant 


que le nombre de nappes est au moins égal à 70 ($ 4.3.1.3 Ÿ ), le nombre de Stanton est donc 
à prendre dans la formule (4.28) avec œ = — 0,7 pour le fluide chaud : 


0,6 
St, = 0,023 | 1 + 6,2 d4 +090] X 14502 | x 500072 pr-07 
1,4 — 0,98 


Le nombre de Prandtl est calculé à la température T, 


me, CE qui donne d’après les 


valeurs de l’énoncé : 
v MCp 310 °xXI212 





Pr.=— 1,10 
a À 0,033 
et l’on obtient : 
St, = 0,024 


Le coefficient d'ëchange h, se déduit de St : 


he = Ste Pe Cpe VË = 0,024 X 0,585 x 1212 x 4,66 


C 


(la vitesse de référence V? étant la même que pour Re, cf. $ 4.3.1.2). 





h.=79,3W/m° K 











Coefficient d'échange he côté fluide froid (dans les tubes) 
On sait que Re p = 25000 ; en T'absence d’information sur la rugosité de la paroi, il est 


raisonnable d’admettre que l’on est en régime « turbulent lisse ». On utilisera donc la formule 
de Dittus et Boelter (4.21) avec œ = — 0,6 (fluide froid), sous réserve que l’on ait L/d > 60, 
ce qui devra être vérifié en fin de calcul : 
si —0,2 p.—0,6 

Stp = 0,023 Re, Pr; 
Pour lair à T,, = 70°C, Pr = 0,7. 
St ; = 0,023 x 25000 7° x 0,7 7° 
St p = 0,00376 


De la nous dëduisons h pë 


i 0 
hys Sts Pf Cf Vf 


A la température T,,; = 70°C, les tables donnent pour lair: pp = 1,02 kg/ m” et 


Cpr = 1015 J/kg K . Ainsi : 
h ş = 0,00376 x 1,02 x 1015 x 10,1 





hş = 39,3 W/m? K 











Coefficient global côté fluide chaud 
Le coefficient global k, côté fluide chaud est ensuite calculé avec (6.6) : 








-1 
E 25 1 
E di Fi 








Er hy 
c f F 
Le rapport des surfaces d’échange est égal au rapport des diamètres D/d = 55/50. 
1 
F , 55 1 
79,3 50 39,3 








k,=24,6W/m°K 








5: 
Puissance 
On peut calculer la puissance par exemple côté froid : 


D = qy (T— Ty )=7,57.10° (120 — 20) 

















D =757 KW 
Efficacité 
D’après (3.12) : 
> p B 757.107 
Qi min (Tee = Tf)  7,57.107 (360 — 20) 


E = 0,295 


NUT 


L’ échangeur est à courants croisés, à un fluide brassë ($ 1.2), avec une seule passe sur 
le fluide brassé. Le débit thermique maximum est celui des fumées, c’est-à-dire du fluide 


brassé. Le tableau 3.1 donne dans ce cas : 


NUT =- Li + 2 Ln(1 - ER) 


UT 97 
9,57 


R 








soit R = 0,79 


dt max 


et finalement : 


NUT = 0,41 


Surface d'échange 


De NUT on déduit la valeur de la surface d’échange 2, côté fluide chaud : 
LEZ - 24,6 X. 
at TO 
E, z 126 m? 





NUT = (0,41 














Faisceau de tubes 


tubes : 


La surface d'ëchange 4. côté fluide chaud correspond au diamètre extérieur des 
C 


D = 55.107? m. Si L, est la longueur totale des tubes, on a donc : 


126 
Z= mDL, soit L, = ner 
Z X 55.107 
L,=729 m 


Appelons n le nombre de tubes nécessaires dans chaque passe pour assurer le débit 


demandé, et S, la section totale des tubes : 





S,=n 
Pour avoir n, il faut calculer préalablement S, à partir de l’expression du débit d’air : 


Im = PJ St Vj 
Dans un calcul de débit, la masse volumique py doit être évaluée à la température 


moyenne de mélange, ici T,ç ($ 4.1.4 et 4.1.5 Y), ce qui donne py = 1,02 kg/m’. En outre, 


V;=V} =10,1m/s,et qmf =7,5 kg/s. Alors: 





e pp 
1,02 x 10,1 
4 

n = dj avec d= 50.107? m 
md 


n = 370 tubes par passe 


soit au total : 





N = 740 tubes 











La longueur de chaque tube (c’est-à-dire la longueur du faisceau) est donnée par : 
L, 729 


= — m 
n 740 














Vérification 
Il y a lieu maintenant de s’assurer que la longueur de chaque tube répond bien à la 
condition L/d > 60 (question 4, calcul de Aș). Or ce n’est pas le cas, puisque nous avons ici 


L/d s 20 : d’après (4.24), Stp doit donc être multiplié par un facteur 7,12. Cela remontera 


de quelques pour cent le coefficient d'ëchange global, et il faudrait normalement procéder à 
une itëration avec cette nouvelle valeur. Mais on est excusable si on renonce à recommencer 
les calculs. 

En ce qui concerne hp, vu le nombre total de tubes, la condition « nombre de nappes 


supérieur à 10 » est réaliste. 
Commentaires 


- La question 2 illustre le genre de raisonnement sur des données incertaines que 
l'ingénieur doit fréquemment pratiquer (ou devrait être apte à pratiquer ....). Bien qu’elles 
soient par essence nuancées, les conclusions auxquelles on aboutit sont appelées à étayer des 
choix. 

- Dans le cas présent, la suite du calcul montre que h, > hp, ce qui a pour effet de 


relever la température de paroi estimée, et donc de faire reculer le risque de condensation 
côté (a) dans la formule « entrées opposées ». Malgré tout, cette solution apparaît moins sûre 
que la disposition « entrées du même côté ». 

- Les caractéristiques thermophysiques des fumées sont fournies dans l’énoncé pour 
raccourcir les calculs. On doit normalement les estimer à partir de la composition des dites 
fumées, avec les tables relatives à O, , No et CO; . 


PROBLÈME N°9: Échangeur tubulaire, faisceau en quinconce 
Énoncé 


On étudie un échangeur air-eau à faisceau de tubes et calandre. L'eau circule dans 
les tubes et lair autour des tubes, perpendiculairement à eux. Les tubes sont disposés 
verticalement, chaque nappe transversale constituant une passe pour la circulation de l’eau. 
La calandre est rectangulaire, de hauteur H „= 2,615 m et de largeur l, à déterminer. 

Le faisceau est en quinconce, de pas triangulaire équilatéral. Le pas transversal entre 
tubes est fixé à Sy = 95 mm. Dans chaque nappe transversale, les tubes sont montés en 
parallèle. Ils sont en acier. Leurs diamètres intérieur/extérieur sont d =31 mm et 
D = 35 mm. Une distance de 0,5 Sy est imposée entre laxe des tubes d'extrëmitë et la 
paroi de la calandre. On se fixe une vitesse débitante maximale pour l’eau (fluide froid) 
Vr=lm/s. 


Les données du problème sont : 


Fluide Eau Air 
Température d'entrëe T fe = 20°C Te = 250°C 
Température de sortie Ty, = 80°C 

Débit Qaf = 12 kg/s Ame = IO kg / s 


1. Calculer la puissance thermique échangée et la température de sortie de l’air. 

2. a) Déterminer le nombre N; de tubes nécessaire dans chaque nappe transversale. 
(préciser quel est dans ce cas le nombre de nappes longitudinales du faisceau de tubes). 

b) En déduire la largeur l, de la calandre et la vitesse y? de lair (fluide chaud) 
dans la calandre vide. 

3. Calculer le coefficient global d'échange ky côté fluide froid. Commenter le 
résultat. 

4. a) Déterminer l'efficacité de l'échangeur, son NUT et la surface d'échange Z ç 
nécessaire (côté fluide froid). 

b) Calculer la longueur totale L des tubes, et le nombre Nr de nappes 
transversales. 

5. La rugosité intérieure des tubes est € = 107? mm. Pour chaque demi-tour du 
fluide (passage d’une nappe à la suivante), le coefficient de perte de charge singulière est 
évalué à Ç =0,4. On comptera Ç, =0,5 pour l'entrée et , = 1 pour la sortie du circuit. 
Calculer la perte de charge totale et la puissance utile P, de la pompe nécessaire pour 
alimenter l'échangeur en eau. 


Solution 


1. 

La puissance de l’appareil se déduit du bilan thermique sur l’eau : 

DSgj(Te Te) 

Avec C,f = 4180 J/kg K environ pour l’eau, le débit thermique unitaire vaut : 
dif = amf Cpf = 12 X 4180 


y = 50160 W/K 


d’où la puissance : 
® = 50160 (80 — 20) W 
PD = 3010 kW 














Le bilan thermique sur l’air s’écrit : 
p = lte (Tee = T )= Ime C pe (Tee cu Tes ) 


En toute rigueur, pour connaître C,.,il nous faut la température moyenne T,,.de lair, 


pe” 
qui est inconnue. Mais les chaleurs massiques varient assez peu avec la tempërature. En se 
basant sur une estimation approximative de 270°C ,ona C pës 1025 JikeK.D'ou: 

Qr = 50 x 1025 

4 = 51100W/K 


On en déduit la température de sortie de lair : 





3 
T =T,- Ê = 250 — 3010.10 
i die 51100 
IE 











La température moyenne T,,. est donc de 220°C , et Tincidence sur la valeur de Cpe 


est négligeable. On en reste là. 


2. 
2a). Nombre de tubes 
Le débit d’eau est donné : 
Les tubes étant en parallèle, leur section totale dans une nappe transversale est : 





2 
Sy = Nr TË = dë 
ST Pl 
Avec Vy-Imis et pp = 10° kg/m*, il faut donc dans chaque nappe : 
N me = 15,9 tubes 


NS KË ER 
zx(31.10 7} x10°x1 
La vitesse de 7 m/s étant un maximum, on doit arrondir à l’entier supérieur : 











Nj 16 tubes 





Le rëajustement qui en résulte sur la vitesse est négligeable, et on peut conserver dans 
la suite des calculs V; =1 m/s. 


Le nombre de tubes d’une nappe transversale est égal au nombre de nappes 
longitudinales, soit ici 16 tubes ($ 4.3.1.1). 


2b). Largeur de la calandre et vitesse de lair 
Une rangée de 16 tubes au pas transversal de 95 mm nécessite un espace de 15 entre- 
axes = 15 Sy , plus 0,5 Sy à gauche et à droite entre les tubes et la paroi de la calandre 


(condition posée par l’énoncé). 


Mais le faisceau étant en quinconce, la seconde rangée est décalée de 0,5 Sy vers la 
droite ou vers la gauche, ce qui fait en tout : 
L=(15+2X0,5+0,5) Sy = 16,5 x 95.107? 
le =1,57 m 














La calandre vide a une section S, = H, l, = 2,615 x 1,57 


S. =4,10 m° 
A la température moyenne T „e= 220 °C (question 1), les tables donnent pour la 


masse volumique p, = 0,715 ket m” , d’où la vitesse débitante dans la calandre vide : 

















yl = mc Z 50 

© PeSe 0,715x4,10 
VË =17 m/s 

3. 


Pour calculer le coefficient d'ëchange global, il faut d’abord connaître les coefficients 
d’échange de chaque côté des tubes. 
3a). Coefficient hç à l’intérieur des tubes 

V, d 
Le nombre de Reynolds dans les tubes est: Re, = —— 
vV 
f 
A la température moyenne de mélange de leau : 


r E 20+80 — soc 
2 





la viscosité cinématique vaut (cf. tables) : V= 0,057.10 7 “mis , d’où : 


1x 31.107? 

AE as EERO 

9057.10 
Re ; = 54400 


Avec ce nombre de Reynolds (< 10 ? ), on doit pouvoir utiliser la formule de Dittus- 
Boelter (4.21), à condition d’avoir un allongement des tubes L/d > 60. Ici la longueur des 
tubes est la hauteur de la calandre, c’est-à-dire H,. On a effectivement : 





Pie os 
d 31 
Le compte est bon, et alors : 
+ -0,2 p,.-0,6 
Sty = 0,023 Re, Pry 


A la température T„f = 50 °C , les tables donnent Pr; = 3,67. 

St ; = 0,023 x 54400 7°? x 3,67 7° = 0,023 x 0,113 x 0,456 

St;=1,185.10 7? 

D’après la définition du nombre de Stanton, avec les valeurs numériques déjà connues: 
h= St; Ps C py Vy T L185.10 °x10°x 4180 x 1 








h;= 4950 W/m° K 








3b). Coefficient h; autour des tubes 


Le faisceau est en quinconce, à pas triangulaire équilatéral. On a donc : 


- pour le paramètre géométrique adimensionnel ($ 4.3.1.3 @ ):  =1/2 
-et pour le pas transversal relatif : er = Br Ka A 
D 35 
er SQ 
Le nombre de Stanton est donné par (4.29b) (dans le cas Y <17 et avec 
a = — 0,7 pour le fluide chaud) : 


E. 0,66 

aTa lie ) Ej RE OP 
er —1 

A la température moyenne de lair 7,,= 220 °C (question 1), Pr, —0,68 et 

Ve = 3,74.10” 3 m? / s. La vitesse de référence pour le nombre de Reynolds est ici la vitesse 


dans la calandre vide ve = 17 m/s (question 2), et la longueur de référence est le diamètre 
extérieur D des tubes d’où : 
VID I7x3510 


Re, à 
Ve 37410 
Re, = 15900 
0,66 
St, = 0,023 l të jeta së ) | 15900 7 0346 x 0,68 “07 
St. = 0,0132 


On en déduit le coefficient d’échange moyen : 
he = Ste Pe Cpe ve = 0,0132 x 0,715 x 1025 x 17 











h. =164,5W/m° K 





3c). Coefficient global k p 

Le coefficient global d’ échange côté froid est donné par (6.6d) : 

2 Z 

e E e EA 

Në Zehe AE hp 
avec ici (§ 6.23%) EplZz d/D et ZylZm” ZdKd + D), et aussi e=(D-d)/2 
= 2 mm. Pour la conductivité thermique de l’acier À = 46 W /m.K . 

E R 2107201 pi 
k 35 164,5 46 31+35 4950 


= 5,38.107? + 0,04.107? + 0,21.107? 


ky 











k; =177,5 W/m° K 








Dans ce total, le terme dominant est le coefficient d’échange côté air. Le coefficient 
d'échange côté eau compte pour moins de 4%. La résistance thermique de la paroi est 
négligeable. 


4. 
4a). Efficacité, NUT, surface d'échange 
D’après la lëre question, q; min = qy - L'efficacité de l'échangeur est donc exprimée 
par (3.11): 
TT 80-20 
Toe- Tj 250 — 20 


ce 


E —0,26 














L'ëchangeur est de type 1-N. Son NUT est donnë dans le tableau 3.1: 


1 2-E(1+R-Ę41+R?) 
NUT = Ln 


TER DCI ER ENTIER") 


Le facteur de déséquilibre vaut : 























R= dt min — 50160 = 0,98 
dimas 31100 
Tous calculs faits : 
NUT = 0,36 
K kri 
Nous savons que NUT = —— = JË 
Qi min Qi min 


La connaissance du NUT permet donc de déterminer la surface d'échange Zp : 





DE NOT a E 0 
ky ea 
Z;=1017 m? 











4b). Nombre de nappes transversales 
Si L est la longueur totale des tubes, la surface d’échange Z f (côté fluide froid, c’est- 
à-dire à l’intérieur des tubes) s’écrit : 
d’où : 
S 101,7 
mx 31.107? 
L = 1045 m 














Dans une nappe transversale, il ya N, tubes, dont la longueur totale est Nj, XH. Le 
nombre de nappes transversales Ny est donc tel que : 
Në L _ 1045 
= N,XH, 16X2,615 
La valeur trouvée est très proche de l’entier 25. On prendra donc, sans avoir besoin de 
faire une correction : 


Nr = 24,97 





Nr = 25 nappes 











5. 
5a). Perte de charge 
Lorsqu’un circuit possède des branches en parallèle, la perte de charge AX est la 
même sur chacune des branches et sur l’ensemble (FEMM, $ 7.6.1). 
Dans une passe, où chaque tube a une longueur H,, d’après FEMM (6.140b) la perte 
2 


V 
<p + et la perte de charge sur chaque singularité s’élève 


de charge en ligne vaut : 4 C; i 





2 
PV} 
2 
De l’entrée à la sortie de l’eau, le circuit comporte Ny passes ; chaque particule fluide 


à: E 


suit un itinéraire qui comporte donc Ny tubes, Ng —1 coudes de demi-tour, une entrée et une 
sortie, avec une vitesse uniforme VF, d’où une perte de charge totale (FEMM. 6.177 et 
6.178) : 


y2 
ax -fse, ND. (ng) E kg ts Jo E 


On a vérifié à la question 3a) que l’allongement H, 7 d des tubes est supérieur à 60. 
Vu le nombre de Reynolds Re ç = 54400, le coefficient de frottement C ; est donné par la 
formule de Colebrook (4.23 et aussi FEMM 7.10) : 


E je 
Le calcul est itératif. Dans la pratique, les valeurs des coefficients de frottement vont à 
peu près de 7 07! à 10 *. Comme point de départ, on peut prendre une valeur moyenne 


C f= 2.107? dans le membre de droite : 


-2 
i — 2,5 Ln)0,285x 1 + ES 


C;/2 31 54400 V2.10”312 





et on obtient comme première estimation : 
GES 


On reprend le calcul précédent en remplaçant la valeur de départ 2.1 07? par cette 
nouvelle valeur, ce qui conduit à une seconde estimation : 


C;=5,20.107* 
Rëinjectons cette nouvelle valeur dans le membre de droite. Le résultat est : 
Ce ROOT 


On est prës de la convergence. Une ultime itëration donne : 





es — 3 
C=5,25.10 











La perte de charge est alors : 
25 x 2,615 


31.107? 


2 
ax =((4x5,2510=* x EKA) OSH To 


3 
—(44,28 + 11,1) — 





AX = 27700 Pa = 0,277 bar 











5b). Pompe 
La puissance utile de la pompe doit être égale à la puissance P, dissipée dans 


l'écoulement (FEMM, $ 7.4.2 et 7.4.4) : 
P,= P, = q, AX 
Cette formule est encore valable pour un circuit avec branches en parallèle, q, étant le 
débit total (FEMM, $ 7.6.1) : 
P, — 12.10” x 27700 
P,=3324W 














Commentaires 


Ce problème prend en compte de nouveaux éléments : échangeur multipasses, 
montage des tubes en parallèle, géométrie du faisceau, perte de charge .... 

- Un calcul approché de la surface d’échange pourrait être conduit à partir du bilan 
d'énergie sur le fluide froid, en appliquant (6.5), soit: B=k;Z;(<T.>-<T;>).Ona 


<T, SET ES 2208 CT ST E 9020 sot E 301010 17177,5 (220.501). 


ce qui donne une surface d'ëchange de 99,75 m” , au lieu de 201,7 m? par la méthode NUT 
(question 4a). L'écart n’est que de 2%. Doit-on en déduire que cette méthode est un 
raffinement inutile ? Certainement pas, d’abord parce que l’écart n’est pas toujours aussi 
modeste (ici, le NUT est petit, cf. Problème 6), ensuite parce que, en accumulant les 
approximations, on finit par élargir sensiblement la marge d’incertitude. 

- Le calcul itératif du coefficient de frottement converge toujours rapidement, même si 
l'estimation de départ est un peu loin du résultat final. 

- On peut être surpris par la faible puissance de la pompe, due en particulier au fait que 
la vitesse de circulation n’est pas très élevée (1 m/s). Mais il faut prendre aussi en compte son 
rendement (rarement supérieur à 0,7) et tout le circuit dans lequel l’échangeur est inséré. Au 
bout du compte, on arriverait facilement à une puissance réelle de 7 à 2 kW. 


PROBLÈME N° 10 : Condenseur 
Énoncé 


On étudie un condenseur de butane, destiné à produire 15 kg/s de butane liquide à 
partir d’une vapeur saturante à la température T,= 50°C et à la pression p = 5,08 bars. 

L'appareil utilisé est un échangeur à faisceau de tubes horizontal, à quatre passes 
côté tubes. Ceux-ci ont un diamètre intérieur d = 15,75 mm et un diamètre extérieur 
D = 19,05 mm. 

On dispose comme fluide froid d’une eau de rivière dont la température maximale est 
de 15 °C en été, et qui ne doit pas subir un ëchauffement supérieur à 10 °C. Elle circule à 
l’intérieur des tubes. 

Les caractéristiques physiques du butane à p = 5,08 bars sont données ci-dessous : 


Température (°C) 20 50 
Chaleur latente de vaporisation Ly (J/kg) 335.105 
Conductivité thermique phase liquide (W/m K) 0,131 0,125 
Viscosité dynamique phase liquide (kg/m s) 0,17.10 01310 
Masse volumique phase liquide ( kg im) 579 542 


1. Calculer le débit de fluide froid q,yç nécessaire dans l'échangeur. 

2. Quel est le débit-masse d'eau q',ç dans chaque tube, le nombre de tubes N 
nécessaire pour chaque passe et le nombre total N, de tubes ? (On se fixe une vitesse 
débitante V = 2 m/s dans les tubes). 

3. Déterminer le coefficient d'échange hy côté fluide froid. 

4. On se base a priori sur une température moyenne de paroi côté fluide chaud 
Te = 30°C. Déterminer quel est alors le coefficient d'échange moyen h, à l'extérieur des 
tubes (raisonner ici comme pour un tube seul). 

5. En tenant compte de la disposition des tubes en nappes verticales (qui diminue h, ) 
et de la résistance thermique de la paroi, on admet finalement un coefficient global d'échange 
côté butane : k.= 730 W/m K. 

Calculer la surface d'échange X nécessaire. En déduire la longueur du faisceau de 
tubes. 


6. La rugosité moyenne à l’intérieur des tubes est € = 0,01 mm. Les coefficients de 
pertes de charge singulières sont donnés : Č çnirës = $ sortie = 1,2 pour chaque tube, compte 


tenu de la forme des têtes de calandre. 
Calculer la perte de charge dans le circuit d’eau et la puissance de la pompe qui doit 
assurer la circulation (en admettant pour celle-ci un rendement global 7 = 0,6 ). 


Solution 


1. Débit d’eau 


La puissance thermique de l’échangeur se calcule en écrivant le bilan sur le fluide 
chaud : 


P = qme Ly = 15 x 335 x 10° 

® = 5000.10? W 

Sachant que l’on doit avoir Th — T to < 10°C, le bilan sur le fluide froid donne, dans 
la condition la plus dëfavorable : 

P Sqy (Ty -Tg )= 10 di 

qy = P/10 = 500.107 W/K = qaf Cpj 

Dans la gamme de température considérée pour l’eau, Cpr = 4180 J/kg K. D'où le 
débit minimal nécessaire : 

Amf = 500.107 14180 





amf = 120 kg/s 











2. Nombre de tubes 


Le débit-masse dans un tube est : 
2 





qt Pj SV py TË V ave V=2ms 


d'n =10°x Ex (15,7510 )?x 2 








G'mf =0,39 kg/s 








Dans chaque passe, les N tubes sont évidemment montés en parallèle, et le débit total 
s'écrit : 

Im EN d'mf 
d’où: N = 120/0,39 
N = 308 tubes par passe 














et avec quatre passes il y a au total : 








N, = 1232 tubes 








3. Coefficient d'échange côté froid 
Pour calculer h on passe par l'intermédiaire du nombre de Stanton, qui s'exprime en 
fonction de Re et Pr. 


Conformément au $ 4.1.5, le nombre de Reynolds dans les tubes sera évalué à la 
température moyenne de mélange de l’eau : 


1 5 o 
Ta PS mi avec T=15 C et Ta = (15 + 10) C 
Tmf = 20 °C 


m, 


D’après (4.13): Reç = a 
Vy 

Les tables donnent la viscosité cinématique de l’eau à 20 °C: 
Vr=0,101.10 7 m°/s 

: 2KTS TRIO 

ef = —5 
O,101x10 
on calculera donc Sr; avec la formule de Dittus - Boelter (4.21) pour le fluide froid : 


St; = 0,023 Re” 92 pr” 05 





= 31200 : l’écoulement est « turbulent lisse » (Re < 10°) | 


Avant d'entamer le calcul, assurons-nous quand même que les conditions de validité 
de la formule sont respectées : 


- nombre de Prandtl : pour l’eau à 20 °C, Pr =7, donc Pr > 0,7. C’est correct. 


- allongement des tubes L/d > 60 7 On ne sait pas. Admettons-le provisoirement. Nous 
verrons à la fin. 


Donc : 

SOINS LION KE 
z -4 

St; = 9.10 





et aussi en revenant à la définition (4.6) : St pa 
PCV 


S'agissant de l’eau : p = 10° kg/m” et Cp = 4180 J /kgK 
hş =9.107fx 107x4180 x 2 








hp =7525 W/m? 








4. Coefficient d'échange côté chaud (butane) 


La condensation se produit à l’extérieur d’un tube horizontal ; le coefficient moyen 
nous est donné par (5.1) : 


EL 0,25 
AE os) ER 
vı D(T,-T,) 
Nous avons comme température moyenne de film : 
_T,+7T, 50+30 
i 2 2 

A cette température, les caractéristiques thermophysiques du butane sont obtenues en 
interpolant les valeurs données dans l’énoncé ; on trouve : 

pi = 554 kg/m? 

À =0,127W/mK 


vı = ÉL —0,2610 5m? / s 


Pi 








= 40C 


En outre : 
Ly 345.107 JI kg 
T,= 50 °C ; T, = 30 °C (valeur admise) 


D =19,05.10 °m 
et le calcul donne : 





h,= 1800 W/m° K 











5. Dimensions de l’échangeur 

La surface d’échange correspondant à la valeur donnée de k est la surface côté fluide 
chaud, c’est-à-dire à l’extérieur des tubes. On a donc, d’après (6.5) et avec le résultat de la 
question 1 : 


D =k, X(T,- Tf ) 730 E, x 30 = 5000.10° 











E, = 228,3 m? 





Z, est la surface extérieure des 1232 tubes, soit en appelant L la longueur d’un tube 


(qui est aussi la longueur du faisceau) : 
228,3 


E, N,TDL etdece fat L= ——_—— 
1232 x x x 19,05.107 





L=3,10 m 











Nous pouvons maintenant vérifier que la formule donnant St a été correctement 
employée dans la question 3 : nous avons ici L/d = 200, donc très supérieur à 60. 


6. Perte de charge 
Dans une passe, les tubes sont montés en parallèle : la perte de charge AX „ est donc 


la même pour un tube et pour la passe (FEMM, $ 7.6.1 et formule 6.178) soit : 
2 
ax, (ec, Emo 
où ZE est la somme des coefficients de pertes de charge singulières. Dans ce problème, 
Z=12+12=2,4: 
Compte tenu de la valeur de Re, le coefficient de frottement est donné par la formule 
de Blasius (FEMM . 7.8a) : 


C = 2 x 0,0395 Re 0” = 0,006 


On peut donc se passer de la rugositë £. C’est une donnée surabondante. 
Pour les quatre passes, la perte de charge totale s’établit à : 


157510 © 
AX = 57000 Pa = 0,57 bar 


2 
ax = ax 006 SHU s 24 ho KE 














D’après FEMM (6.167a) et (7.73b), la puissance dissipée dans le circuit est : 
P, =q, AX 

En admettant que l’entrée et la sortie du circuit sont au même niveau, la puissance 
utile fournie par la pompe doit être égale à la puissance dissipée dans le circuit. On a donc à 


partir de FEMM ($ 7.4.2, $ 7.5.2.1, et relation 7.55): 

















P AX m 
per dh avec q, = MË =0,12m°/s et  m=0.6 
n n P 
P = 11400 W 
Commentaires 


1. L'ënoncë contient certaines informations qui ont pour but de faciliter la résolution 
du problème, mais qui doivent normalement être estimées par l'ingénieur : il s’agit du 
coefficient moyen d'échange k, et de la température moyenne de paroi T, . 

En fait, ces deux grandeurs sont interdépendantes puisque la formule qui donne k, 
contient T,. Là comme dans d’autres exemples cités, c’est une procédure itérative qui est à 
recommander : 

- première estimation de T, 

- calcul de h}, k, et du nombre de tubes 

- choix de la disposition du faisceau et en particulier du nombre de tubes dans une nappe 
verticale, et réévaluation de k,. par (5.2) 

- nouvelle évaluation de T, et réitération du calcul. Souvent, la convergence est obtenue dès 


cette seconde étape. 
2. Pour le calcul de C fe la formule de Colebrook (FEMM, 7.10) est plus sûre, mais 


Blasius est plus rapide et donne ici un résultat très voisin. 
3. La valeur de h relative à la condensation du butane est nettement plus faible que ce 
que l’on obtiendrait avec de l’eau. 


PROBLÈME N° 11: Échangeur compact 
Énoncé 


On aborde l'étude d'avant-projet d’un échangeur compact à courants croisés, avec un 
fluide brassé. 

L'ëchangeur est constitué d’une batterie de tubes parallèles, reliés par des ailettes 
perpendiculaires aux tubes et jouant le rôle d'entretoises (fig. 1: vue de face, et fig. 2: coupe 
longitudinale d’un tube). On admettra que les tubes sont de section rectangulaire (fig. 3 : 
coupe transversale d’un tube), disposés en quinconce (fig. 4 : vue de dessus). 

II s’agit d’un échangeur air-eau. L'eau circule dans les tubes, l’air autour des tubes, 
entre les ailettes. Les dimensions du caisson sont : 1 m (largeur) X I m (hauteur) X 0,5 m 
(profondeur). Les autres dimensions sont portées sur les figures, en mm. L'appareil comporte 
1786 tubes et 500 ailettes. Les caractéristiques aux entrées sont les suivantes : 


- pour l'eau: Ty, = 12°C; qyp = 28 kg/s 
- pour l'air: T = 125°C ; qme = 14 kg/s 








Fig. P11.1 Fig. P11.2 
TË 
—— air 
: © 
— AT 
1 E) 
— air 
| ailette 
Dë Pe ët vën mëmë ën ad ËS J 
0,25 © 
0,5 m 


Fig. P11.3 Fig. P11.4 


1. Déterminer le coefficient d'échange hy côté eau. 

2; 
2.a) Calculer le diamètre hydraulique pour l'écoulement de l’air (la section à prendre en 
compte est la section frontale de passage de l’air, voir fig. 1). 
2.b) La section totale de passage de l’air est S, = 0,782 m°. Evaluer le nombre de Reynolds 
Re, côté air. Quelle est l’influence de la pression de l’air ? 


2.c) L'abaque de la figure 5 donne St Pr?” 


en fonction de Re,. La température de 
référence est la température moyenne de mélange T,... Calculer le coefficient d'échange h, 


côté air. Comment est-il influencé par la pression de l’écoulement ? 

































































Fig. P11.5 


3. La surface d'échange X, côté fluide chaud est de 489 m”, dont 418 m° pour les 
ailettes. On admet pour ces dernières une efficacité € = 0,86. Déterminer la conductance 
globale K de l’échangeur (on néglige la résistance thermique des parois). 

4. Calculer l'efficacité de l échangeur 

5. Calculer les températures de sortie pour l'air et l’eau. Vérifier si les températures 
moyennes adoptées en cours de calcul sont suffisamment proches des températures calculées. 


Solution 


1. Coefficient d'échange côté eau 

La première chose à faire est d’estimer approximativement les températures de 
référence pour évaluer ensuite les caractéristiques thermophysiques des fluides. S'il s’avère à 
la fin du problème que cette estimation est trop imparfaite, on recommencera en se basant 
cette fois sur les résultats du premier calcul. 

Ecrivons déjà le bilan enthalpique global : 


mc C pe (Toe g Tes )= mf Cof (Ty = T to ) 


Les chaleurs massiques varient peu avec T. Prenons C,,.= 1007 J/kg K pour lair et 
Cf = 4180 J/kg K pour l’eau, soit compte tenu des valeurs imposées aux débits : 
14 x 1007 x AT, = 28 x 4180 x AT; 





AT, 
ce qui donne —— = 8 
AT 
Un ordre de grandeur plausible serait AT, = 90°C et AT, = 11°C, ce qui respecte à 


peu près le rapport de 8. Cela correspondrait à des températures de sortie T., = 35°C et 


T = 23°C, d’où des températures moyennes de mélange ($ 4.1.5V): 


E o TI 23 


t së = 80°C; Th E = 17,5 °C 





Nous partirons sur ces bases-la. 


La section du tube est rectangulaire ($ 4.2.1.3, 2b = 2 mm, l = 18 mm) et Re est 
donné par (4.17) : 
V; Dyf N 2x2bl 2x2xl& 
tj ou hf — — 
It 2b 18 + 2 


= 3,6 mm 





Re; = v, 
La section intérieure d’un tube est: S = 18 x 2 = 36 mm”. Elle est parcourue par un 
débit : 
ms ae. : 2e 
Lg nb. de tubes 1786 
d’où une vitesse : 
__15,68.10 7 
= 310 
Re ; est calculé à la température moyenne de mélange. A Ty = 17,5°C, les tables 


= 15,68.10 7° m? / s 


= 0,435 m/s 


numériques donnent pour la viscosité de l’eau : 
Vr=1,06.10 mis 
3,6.10 7° x 0,435 
1,06.10 75 





On en déduit : Re = 


Reç = 1480 


L'ëcoulement est donc laminaire. 


La formule à utiliser pour calculer St est en principe (4.18) ou (4.19). Pour trancher, il 
faut comparer L7 Dyç et 0,014 Pe. D’après (4.17) : 





Pe; = 
f 
a 
avec dy = 0,142.10 7° mis pour l’eau à Tmf = 17,5°C , ce qui donne : 
Pa 3 
Pe ç = 11.10 


Ici, LzIm,soit: 


L 1000 
Dpp Pe 3,6 x 11000 


C’est donc la formule «tubes longs » (4.19) qui s'applique, sous réserve que les 
conditions sur Re, Pret 4/ 4, soient bien vérifiées, ce qui est le cas : 


8 23 0,14 
sr, = #23 EJ 
Pes \ Hp 


En ce qui concerne / U, , on doit s'attendre à une température moyenne de paroi T, 


= 0,025 > 0,014 


assez proche de T,,; puisque le coefficient d'échange côté eau sera en principe beaucoup plus 


élevé que côté air. En admettant T,=T,,; + 10°C, on trouve 4/1, = 1,07/0,88 d’où une 


0,14 
) 


correction ( 4 / = 1,03 que l’on pourrait à la rigueur négliger. 
H/ Hp q P 8 818 


On obtient finalement : 
8,23 


re 
f 11109 





X 1,03 


Stp = 0, 77.1073, et alors : 
hp = St pC, V =0,77.107 °x 10° x 4180 x 0,435 
hş = 1400 W/m° K 














2. 
2.a) Le nombre de Reynolds se calcule dans un canal élémentaire, dont la section est 
11,5 x 1,9 mm (fig. 1). Le diamètre hydraulique est donné par (4.17) : 


e DAKO QA KOKA 




















NT 142b (1154 1,9)1055 
Dy, = 3,26 mm 
V.D 
2.b) Nous avons ici: Re, = — he et y, =% 
E C 


En admettant une bonne équirépartition des débits dans les canaux, la vitesse est la 
même partout ; elle peut être calculée directement à partir du débit et de la section totale : 


V; — fme , soit: 
Pe Se 
Re, = ane Prc 
He Sc 


On ne nous précise pas la pression moyenne de l’air, mais cela est sans importance car 
la viscosité dynamique dépend très peu de la pression (propriété indiquée dans les tables 
numériques). Comme c’est le dëbit-masse qc qui est fixé, Re, est donc indépendant de la 
pression. 


2.c) Le calcul doit être fait à la température moyenne de mélange, soit T,,.z 80°C . Alors : 


U, 2,08.10 7° kg/m.s 





14 x 3,26.10 7° 
Re, = TES T 
2,08.10 7 © x 0,782 
Re, = 2800 











L'abaque donne approximativement : 

St, Pr’? — 6,510” 

Pour l’air, le nombre de Prandtl est peu dépendant de T : Pr = 0,7 
St.=6,5.10 °x0,7 7?" 











St. = 8,24.10 
On a aussi : 
ho q 

St, 5 —— avec V.=- = G, vitesse massique (4.10 

E pe Ca V, Pe C Sa C q ( ) 
De ce fait nous n'avons pas besoin de p, et le coefficient h, ne dépend pas de la 

pression : 
he = Ste Cpe EH avec C,.= 1007 J/kg (question 1) 
c 

Adoncque : 
pa 8,24.10 7° x 1007 x 14 

i 0,782 











h.=149W/m° K 





3. 
Les ailettes sont situées côté fluide chaud. La conductance globale est donnée par (6.8) 
et (6.6a) : 





-1 
1 1 
= + 
h (ZL te2,) Zyhy 
avec X, = 418 m° ; Z; = Z, — X, (surface latérale des tubes) = 71 m’; € =0,86. 
Côté eau, la surface d’échange est la surface intérieure des 7786 tubes : 
Z'} = périmètre XI mx 1786 =2(18+2) 10” x 1786 


Priim 
Revenons à K : 





-1 
K = 1 A 1 
149 (71 + 0,86 X418) 71,45 x 1400 





K = 39370 W/K 











4. 
En revenant à la question 1, on trouve pour les débits thermiques unitaires : 


Gr = 14 x 1007 = 14100 W/K 


dy = 28 x 4180 = 117.10°W / K 
Donc, q; min = Qe = 14100 W / K 
L'air chaud qui circule autour des tubes est brassé. Si l’on se reporte au tableau 3.1, on 
est à la rubrique « échangeur à courants croisés, un fluide brassé à q; min » : 
E=1-exp(-T/R) avec 1°=1-exp(- R NUT) 
_ tmin _ 14100 


R — 
dimax 117000 





— 0,12 


NUT est alors calculable par (3.14b) : 




















e 070 

dimin 14100 
NUT - 2,79 
T = 1— exp(— 0,12 x 2.79 ) = 0,285 
E —1-exp( 0285 

0,12 

E = 0,907 
5. 


5.a) Température de sortie d'air 
Puisque q= q; min» l efficacité s'exprime selon (3.10) : 
E= Tee në Tes 
Tee — T to 
On en tire : 
TËS IT SET, T  YSAZI — 0,907 (125 = 12) 
Tey = 22,5 °C 














5.b) Température de sortie d’eau : 

Le bilan enthalpique nous dit que : 

dtc (Tee QE Tes )= dif (Ts z T to ) 

PS IZa I (125 —22,5) 
117000 

T = 24,5 °C 

















9.c) Vérifications 
- Pour la température moyenne de l’eau, nous avions choisi 77,5 °C. Avec les valeurs 
calculées elle est de (12 + 24,5 )/2 = 18,2°C. La différence est insignifiante. 


- Pour la température moyenne de l’air, l’estimation avait donné 80 °C. Nous sommes 
à (125 + 22,5)/2 = 74°C. L'écart n’a pas de répercussion significative sur les constantes 
thermophysiques. 

- Quant à l’estimation de T, (question 1), son ordre de grandeur est compatible avec 


les valeurs calculées de h, et hy. Sa révision ne changerait pas grand-chose à la correction en 


MI Ulp qui reste faible de toute façon. 


Commentaires 


Le problème s’inspire d’une application présentée dans l’ouvrage de J.F. Sacadura 
[Initiation aux transferts thermiques]. 

L'appareil est assez performant. Son efficacité est élevée, et on observe que la 
température de sortie du fluide froid (24,5 °C) est supérieure à la température de sortie du 
fluide chaud (22,5 °C), situation très favorable (cette éventualité a été signalée dans le cas des 
échangeurs à contre-courant, $ 2.3.5). En pratique, cet excellent résultat est probablement à 
nuancer un peu, du fait notamment que les débits ne sont jamais parfaitement équirépartis 
(solution question 2.b) dans les tubes ou dans les différents chemins de passage de l’air. 

Concernant la section rectangulaire des tubes, on notera que malgré un aplatissement 
marqué (2/18 mm) l’approximation 2 b KI est à éviter : elle donnerait un diamètre 


hydraulique de 4b = 4 mm alors que l’on trouve Dy, = 3,6 mm avec la formule complète 


(question 1). Même remarque pour la section de passage de l’air (question 2a). On peut 
signaler également que, en réalité, les extrémités des tubes seront arrondies pour améliorer 
l’écoulement et diminuer les pertes de charge côté air. 


PROBLÈME N° 12 : fb sanë 


Énoncé 


On dispose de deux échangeurs identiques et de deux sources d’eau (une source 
chaude et une source froide) dont les températures Tee, Ty et les débits qe, 4mp sont 
donnés. 

L'objet du problème est de choisir, entre plusieurs dispositions des échangeurs, celle 
qui permet d'obtenir la valeur de Ty, la plus élevée. 


1. Avec un échangeur seul, montrer que quelque soit son modèle, si qy min = 4w ona: 
Tes — T te Z 1- E 
Te-Tgy 1-RE 


2. Les deux échangeurs sont placés en série comme l'indique la figure P12.1. 
2.1. En raisonnant sans faire de calculs, montrer que les deux appareils ont la même 
efficacité (désignée par E’). On admet que l'influence de la température sur les constantes 
thermophysiques est négligeable. 


2.2. Soit E, l'efficacité de l’échangeur équivalent à l’assemblage réalisé. En admettant 
Qt min = dre» et en utilisant le résultat de la question 1, exprimer E, en fonction de R et de E”. 


Déterminer Ty, en fonction de Toe, Tf, E, et R. 





Fig. P12.1 


3. Les ëchangeurs sont maintenant disposës en sërie sur le fluide chaud et en 
parallèle sur le fluide froid (fig. P12.2). On ne précise pas ici quel est le q, minimum. 


3.1. Montrer sans calculs que les deux appareils ont la même efficacité ( E', côté chaud ou 
E' ç côté froid). 

3.2. En raisonnant sur le circuit série, c’est-à-dire côté fluide chaud, exprimer T,, puis Ty, 
en fonction de E', et des données. 


4. Ceci n’est pas une question. 
En permutant les circuits (fig. P12.3) on obtient : 


' ' ' 2 
Th=Tee E (2 E ç )tTy(l- Ep) 





Fig. P12.2 


5. Les deux appareils sont des échangeurs à courants croisés, avec fluide chaud 
brassé à une seule passe, de surface d'échange X = 170 m? . On a déterminé 
expërimentalement le coefficient global d'échange k. Pour qmo 74 kg/s et qy <10 kg/s : 

k = 52 (qy + 0,8 Qme) +65 en W/m°K 

On donne : Toe = 90°C; Ty, - 10°C 

Déterminer l'efficacité d’un échangeur seul dans les trois cas suivants : 

a) qm = 9,6 kg8/ 5, Qme = 8,2 kg/s 


b) qm = 4,8 kgl si me = 8,2 kg/s 
C) Amp = 9,0 kal ss qe =4,1kg7s 


Calculer Ty, dans le cas a). 





Fig. P12.3 


6. On rappelle les deux sources d’eau dont on dispose : 
- eau chaude : T. = 90°C; qme = 8,2 kg/s 
- eau froide : Ty, 5 10°C; qy = 9,6 kg/s 

Calculer Ty, dans les trois dispositions suivantes (on désignera avec le symbole ' les 
caractéristiques relatives à un seul échangeur) : 


I - en mettant les ëchangeurs en série 
II - en les mettant en série sur l’eau chaude et en parallèle sur l’eau froide (le fluide brassé 


est le fluide chaud) 
III - en les mettant en série sur l’eau froide et en parallèle sur l’eau chaude (le fluide brassé 


étant toujours le fluide chaud). 
Conclure 


Solution 


1. 
Sachant que q; min = 4x » d'après (3.10) l’efficacité et le facteur de déséquilibre 


s’écrivent : 


E — Te” Ts pd që 
Lite Qy Tears 
On vérifie alors aisément que : 
ps Te Tes 
1- E le T to Tej” T to 


= = CQFD 
I-RE 1 Ty E Te Tee Tas To= Ti 


Tee — Tes Tee — T to 








2. 
2.1) Les deux appareils sont identiques : ils ont en particulier la même surface d'ëchange Z. 


Etant traversés par les mêmes débits chauds et froids, ils ont aussi le même coefficient global 
d’échange k, et par conséquent le même R et le même NUT. Leurs efficacités sont donc 


égales. 


2.2) D’après la question 1, avec les notations de la figure 1 : 
Ty Tea 1-E' 
Te-Tp 1-RE 

Tes Te _ 1-E 

Ty -Tp 1-RE 

Effectuons le produit de ces équations membre à membre : 

PART (RER 

TT (I-RE) 

Si l’on considère l’échangeur équivalent à l’ensemble, pour lequel q; min est toujours 








égal à g,, On a aussi : 


T, 


o Lh të l-E, 
ete IKE 


S 
Des deux relations précédentes on tire : 








E (1- RE el 
“I= RE) -ROEF 











Puisque q; min = 4tc dans le réseau, l’efficacité de l’ échangeur équivalent s'ëcrit : 


Tente 
E, = ce CS d’où : 
Tee a T to 


Tes = Tee ~ Es (Tee 7 T po ) 


Nous disposons également du bilan global du réseau : 
di (Ts a T to ) = qr (Tee ~ Tes) 
d’où l’on peut extraire Ty, : 


q 
Th — Tet KE (Te Tes ) 


dif 


En remplaçant T, par son expression ci-dessus, on obtient la température de sortie du 


fluide froid : 











Tp =T p+ RE, (Toe Tp) 





3. 


3.1) Le raisonnement est analogue à celui de la question 2 : les échangeurs sont identiques et 
travaillent dans les mêmes conditions (débits, coefficients d’échange etc). Par conséquent ils 


ont le même R et le même NUT, donc la même efficacité. 


3.2) L'efficacité côté fluide chaud s’écrit, selon l’ échangeur considéré : 
ye Tee — Ta Z Tal D Tes 
a 
Tee- Tie La T pe 


ce e 





Ces deux équations constituent un système à deux inconnues 7, et T,., 
l’une d'elles : 


Pam ET 
Te] = E”, (Ta E T to 14 Tes soit Te 3 E 
te 
et reportons dans l’autre : 
e Ee Tg 
Ecl To= Lfe) = 7 cl he —— 
L= Ey 


On obtient ainsi la température de sortie chaude : 











Tes = Le EL = ke ) (Tee = T i ) 





. Tirons T,; de 


Le bilan enthalpique global du réseau donne ensuite la température de sortie froide : 
di (Ty = Te ) = qre (Tee ~ Tes ) 
En remplaçant 7,. il vient : 





dte 
di 





Ty= Ty + ME E (2 — E, )(Tee-Th) 











4. 
La question n’est pas posée mais le calcul serait analogue en raisonnant encore sur le 
circuit série, c’est-à-dire cette fois-ci sur le circuit froid. 


5. 
5.a) La chaleur massique de l’eau dépend très peu de la température ; de ce fait les C, sont 


pratiquement identiques sur les deux circuits, et l’on voit que : 
dt min — dr du fluide brassé = die 


Il y a une seule passe sur le fluide brassé. Du tableau 3.1 on tire alors : 
E =1-exp(-T/R) avec  Z'=1—-exp(- R NUT) 

et toujours : 
NUT —KkZI qy min 


Les donnëes permettent de calculer k : 
k — 52 x(9,6 +0,8 x 8,2)” +65 


k— 274 Mim K 
Dans la plage de température indiquée on prendra pour l’eau : 
Ch = 4190 J/kg K 


Qi min = Ame Cp = 8,2 X 4190 = 34360 W/K 





On en déduit NUT et R : 

NUT = 274 x 170 = 1,35 
34360 

R = dte - dme = 8,2 = 0,854 


If Im 206 
De la on tire : 
I = 1 — exp(— 0,854 x 1,35 ) = 0,684 
E = 1 — exp( — 0,684/0,854) 





E = 0,55 











En se référant à (3.10) on obtient enfin : 





Ta = 47,6 °C 











5.b) Maintenant, q; min = q; du fluide non brassé = qy . Le tableau 3.1 nous dit que : 


1 

E=-Ù-exp(-RT)] avec l'=1-exp(- NUT ) 
Nous avons donc cette fois : 
k=52(4,8+0,8 x 8,2 )°° +65 = 240 W/m° K 

qi min = Amf Cp = 4,8 X 4190 = 20110 W / K 

kX 240x170 








NUT = = 2,03 
qini 20110 
R-— MË 4S 0,585 
mc 8,2 


ce qui donne : 
IT = 1 — exp(- 2,03 ) = 0,869 


1 
E = —— |] — — 0,585 x 0,869 
F [1 — exp( )] 








E = 0,68 








5.c) La situation est la même que dans 5.a) : Gy min = q; du fluide brassé = qy. 
La succession des calculs est analogue : 
k =52 (9,6 +0,8 x 4,1)” +65 =25]1,5 W/m? K 





nons UU e 
4,1 x 4190 
RER (as 
9,6 
I = 1 — exp(— 0,43 x 2,5 ) = 0,658 


et enfin : 
E = 1 — exp( — 0,658 / 0,43 ) = 0,78 








E = 0,78 








6. 
6.1. Montage en série 

Avec les valeurs numériques données, cette disposition correspond au cas 5.a). On a 
donc pour chaque échangeur : 

E'=0,55 ; R' = 0,854 

D’après la question 2, l'efficacité de l’ échangeur équivalent a pour valeur : 
(1—0,55 x 0,854) —(1-0,55 )£ 


E= — L 
(1 -0,55 X 0,854)? — 0,854 (1 — 0,55) 


E, = 0,73 
et toujours d’après 2. : 
T s = 10 + 0,854 x 0,73 x 80 


Ty, = 60°C 














G.I1. Série sur circuit chaud et parallèle sur circuit froid 
Les données correspondent cette fois au cas 5.b). 


Ame = 8,2 kg/s 
Le débit froid dans chaque échangeur est égal à la moitié du débit total : 
; 1 9,6 
Tmf = 5 Inf Ro ES 
et par conséquent : 
gr min — di 
D’après 5.b) on a donc : 
R' =0,585 
E' = 0,68 


A partir de la question 3, on a : 
Tp= Tp +E E ((2-E', )(Tee= Tp) 
If 

où qre et qy sont les débits globaux dans le réseau. 

Pour calculer T, nous avons besoin de E”, (efficacité d’un échangeur côté circuit 
série, c’est-à-dire circuit chaud). Mais ici E' = E” ç puisque q; min = Ty et: 

ET R'E'ç = 0,585 x 0,68 = 0,4 
donc : 





82C, 
Ty = 10 + x0,4X(2—0,4)x 80 
9,6 C, 














6.III. Série sur fluide froid et parallèle sur fluide chaud 
La situation est celle du cas 5.c). Pour chaque échangeur : 


1 
T'mf = 9,6 kg/s de E no 
E' —0,78 
La formule de la « question 4 » qui donne Ty, fait intervenir E'ç : 
' $ ' 2 
Or, ici E' = E', puisque q; min = g'- Donc E, S RE, : 
[1C 
Le —— — 0,427 
dy 9,6 Ch 
E'ç = 0,427 x 0,78 = 0,333 


Re Le = 


et finalement : 
Tg = 90 x 0,333 x (2 — 0,333) + 10 x (1 — 0,333) 
Ty 54,3 °C 














Conclusion : 
- avec un échangeur seul : Ty, = 47,6 °C 


- avec deux échangeurs en série : Ty, = 60 °C 
- les deux montages en série-parallèle sont presque équivalents : Ty, = 53,7 ou 54,3 °C. 


Le montage en série est ici nettement meilleur. 


Commentaires 


Le problème reprend sous une forme plus simple, avec deux échangeurs seulement, les 
raisonnements du chapitre 7 ($ 7.2 et 7.3). Mais on ne doit pas, à partir des résultats obtenus, 
déduire une règle générale sur la supériorité de tel ou tel assemblage. 

- Le fait de donner une loi phënomënologique k = f(g) permet de calculer k 
rapidement et donc d’abréger la solution. 

- Les questions sans calcul (2.1. et 3.1.) déroutent souvent les étudiants, mais elles 
permettent de voir si les notions essentielles ont été comprises. 

- Dans les questions 5.2. et 5.3., l’erreur à ne pas faire est de confondre le débit total et 
les débits partiels sur les branches en parallèle, ou de confondre E”, et E'ç. 


- On doit bien noter que sur le circuit en parallèle, les températures de sortie des deux 
échangeurs ne sont pas les mêmes (7,, + T,9, $ 7.3). Pour calculer les efficacitës en revenant 


à la définition Ø/ Pmax il faudrait connaître ces températures. C’est pour cela qu’en série- 
parallèle on raisonne sur le circuit série avec les efficacités relatives. 


PROBLÈME N° 13 : Sinë 
Énoncé 
On étudie un échangeur tubulaire en épingle, selon le schéma ci-dessous (un seul tube 


a été représenté en pointillés). Ses caractéristiques sont E, R et NUT. Le fluide chaud circule 
à l’intérieur des tubes. 





Fig. P13.1 


Ce modèle d'ëchangeur peut être considéré comme équivalent à l'assemblage de deux 
demi-échangeurs de mêmes dimensions, l’un (N° 1) fonctionnant à co-courant, l’autre (N° 2) 
à contre-courant (fig. P13.2). 


Te Thos 












FIG. P13.2 


Les données du problème sont : 
; ; os SË 1 
- les débits thermiques unitaires (avec la condition q, < E qy ) 


- les températures d'entrée T,, et Ty, 


- les diamètres d, D, et la longueur totale L de chaque tube (L, = L, = L/ 2). On admettra 


que les tubes sont tous identiques, et on négligera les problèmes liés aux coudes qui assurent 
le demi-tour. La calandre est isolée à sa surface externe. 


I. Justifier la relation suivante: Ty, S (Ty +Ty2)/2. Exprimer R; et R, en 
fonction de R , ainsi que NUT, et NUT, en fonction de NUT. 
2. Calculer E, et E, en fonction de R et de NUT. Calculer ensuite T., (température 


de sortie chaude du demi-échangeur N°1) et T., puis l'efficacité E de l'échangeur complet 


en fonction de E,, E, et R (comparer avec le résultat de la théorie générale des réseaux, 
ch.7). En déduire Ty,. 
3. L'appareil est un échangeur eau-eau. On donne : 
Tie= 90°C ; TR=20°C ; qy=82m/h; q%e=19,80m /h 
faisceau de 60 tubes en parallèle; L=5m ; d=15 mm; D =18 mm 
pas carré Sy = S; =28 mm ; diamètre de la calandre D,.= 0,40 m 
Effectuer sans calculs une première évaluation des températures moyennes de 
mélange T pe et Tr. 
Calculer le diamètre hydraulique D, et le nombre de Reynolds Re; de 
l'écoulement en calandre. Déterminer le coefficient d'échange hy autour des tubes. 
4. Calculer le coefficient d'échange h, à l’intérieur des tubes. Déterminer le 


coefficient global d'échange k (on raisonnera comme si la paroi était plane, et on négligera 
sa résistance thermique). 
5. Calculer le NUT et l'efficacité E, puis les températures de sortie Ty, et T,,. 


Solution 


1. 
la). Le débit froid dans chaque demi-échangeur est la moitié du débit total : 


1 
Imf1 — Amf2 5 2 Tmf 
Donc la température de sortie Ty, est la moyenne des deux températures de sortie T g7 


et Tsz š 


; 1 
1b). Dans le demi-échangeur N°1, ona: qp; = 2 q À Gjer = de 


I : 
Mais on sait aussi que q, < 5 qıy - Donc Gjej < q1- En conséquence : 











dtl min 7 tcl — ic 





R, = Ti1 min. del 
dtl max dfl 
R; = 2 dte 
dif 
Pour l’échangeur complet : q; min = 4i et R = Ge. 
dif 














Le NUT du demi-échangeur N°1 est : 
kı Z; 


dtl min 


NUT, = 


Le coefficient d'échange est indépendant du sens de circulation des fluides, donc 
kj = k, = k coefficient d'ëchange global de l’ échangeur complet. 


La surface d'ëchange X, est la moitié de la surface totale Z. 





Pa 
NUT, = KZ 
2 Ge 
Z pa 
Le NUT de l’échangeur complet vaut : NUT = re 
dt min dte 
D'où 





NUT, = Z NUT 

















Le raisonnement est identique pour le demi-échangeur N°2 : q,2 min = 4te» €t : 














; NUT, = = NUT 








2. 
Za). Efficacitës 


E; = L 





Le demi-échangeur N°1 fonctionne en co-courant. D’après le tableau 3.1 : 


TËR {1 — exp[-(1+ R; ) NUT, |} 


1 1+2R 
E; = 1 — exp| - —— NUT 
1+2R 2 


Le demi-échangeur N°2 est à contre-courant. Toujours d’après le tableau 3.1, son 
efficacité est : 

















+ 1- exp[-(1- R, ) NUT, | 
“T= R exp (i= R; NUT] 


1-2R 
1 - ep - nur | 
2 
E, = 


1-2R evj- l E ur | 























2b). Températures de sorties chaudes 


Dans le demi-échangeur N°1, Qyj min = dre - Alors, d’après (3.10) : 


’ ce c2 
Tee — Te 
To = T. 


ce — E; (Tee z Te) 


Dans le numëro 2, c’est pareil : q,2 min = Gj. Mais ici, la température d’entrée chaude 
est T,9. De ce fait : 


To — Tes 
To — T to 
Tes = To z E, (Ka i Te ) 
E Tee a El a T to De E, (fr. zi Er (Tee H T te ) J T te 





T, 


CS 


ST, 


ce 








See) CB -EE + Ez) 





Z 92 


Dans l’échangeur complet, q; min = 4, et l efficacité s’écrit : 
z Tee — Tes 

C 

Tee = T to 


E —Ej —- Ej E, tË, 


E-E 














Les deux demi-ëchangeurs sont assemblés en sërie-parallële, le circuit série étant celui 
du fluide chaud. Le résultat est conforme aux expressions (7.29) $ 7.3.1, ou (7.32a) $ 7.3.3. 


On a également, toujours d’après (3.10) : 


d’où l’on tire : 
TieTer RE(Toe—T fe) 














3a). AT moyennes 

Vu le rapport des débits, qui est voisin de 4, le AT; de l’eau froide (4qmax) vaudra 
environ 4 du AT, de l’eau chaude. On peut admettre, pour commencer le calcul AT; = 10°C 
et AT, = 40°C , ce qui donnerait pour les moyennes : 

Le 29 Ci: TË TO 


m, 


3b). Diamètre hydraulique de l'écoulement en calandre 
Le faisceau est longitudinal ($ 4.3.2) ; D, est donné par (4.31c) : 








Di=4 Lsro — 18.1073 = 0,0554 — 0,018 
ZX18.10 7 
D, =0,0374 m 











3c). Nombre de Reynolds 
Vs D} 
“ai 





Rep = avec V, = vitesse débitante du fluide froid, qu’il faut calculer. 


La section de passage du fluide froid est : 

S = section calandre — section tubes 

Le faisceau comporte 60 tubes, mais comme ils sont en U il y a 120 tubes dans une 
section droite de la calandre : 








D2 2 

Sp = KS — 120 ZPL =Z (0,40? — 120 x 0,018? ) 
4 4 

S p = 0,095 m° 


La vitesse débitante est définie par : 

f z= —— = ——— K— 

Ss 3600 0,095 

Vs =0,24 m/s 
La viscosité de l’eau à T,,; = 25°C est v p = 0,92.10 “mis. 

0, 24 x 0,0374 
Re} = 5 

0,92.10” 

















3d). Coefficient d'ëchange cotë froid 
Pour le nombre de Stanton externe des ëchangeurs tubulaires a courants parallëles, on 


dispose de la formule (4.32), valable si 5.1 0 < Re < 10° , Soit avec le fluide froid : 
2 -0,18 p..—0,6 
SË p — 0,026 Re, Pr; 
A 25 °C, le nombre de Prandtl de l’eau est: Prş =6,4. 
St ; = 0,026 x 9760 7°" x 6,47” 
Në -3 
Stp = 1,63.10 
On calcule le coefficient d'ëchange Ap en revenant à la définition de St : 
hp= Sts Pj Co Vp 
et en prenant C „p = 4180 J/kg K pour l’eau : 
hş = 1,63.107 7x 10° x 4180 x 0,24 











hp =1635W/m° K 





4. 
Aa). Coefficient d'ëchange interne 
Dans chacun des 60 tubes, le débit est : 
ve, 19,80 


PO “60 60x 3600 


4,0 = 0,0917.10 7° m°/s 
d’où une vitesse de circulation : 





= 40 __4X0,0917.107* 

er 2 di 2 

zd”14  nx(i5.105) 
V,—0,52 m/s 





Le nombre de Reynolds vaut : 
Vd 
NG 
On a admis sur le circuit d’eau chaude une température moyenne T,,. = 70°C , d’où la 


Re 


viscosité V, = 0,42.10 7 m’ /s. 





0,52 x 15.107? 
Re, = 5 
0,42.10” 
Re, — 18600 


La formule de Dittus-Boelter (4.21) pour le fluide chaud s'applique si L/d > 60, ce 
qui est largement vërifië (5000715) : 

St, = 0,023 Re; °°? Pr: 7 

A Tae =70°C , Pr, = 2,62 

St, = 0,023 x 18600” 0“ x 2,62707 

St, = 1,64.10 7? 





Le coefficient d'ëchange est : 

h,= Ste Pe Cpe Ve 

La chaleur massique varie peu avec T. On gardera celle du circuit froid. 
h, = 1,64.10 7° x 10° x 4180 x 0,52 


h.= 3565 W/m° K 














Ab). Coefficient d'échange global 
La formule 6.2a (paroi plane) s’écrit en l’absence de résistance : 


Report a = (0,28 + 0,612)10 7° 
k hë hy 3565 1635 


C 








k=1120W/m° K 











5. 
5a). NUT 

Pour l'échangeur complet, puisqu'on raisonne comme avec une paroi plane, on a 
d’après (6.10 a et b) une conductance globale : 

KOSH, 
avec 


D + d 
Zii =m 2 L çotale 





Alors : 


NUT = K _ KZ» 


dt min dte 








Pour les 60 tubes : 
ra= ax x 1077 Je 5 «60 — 15.56 m? 


Il faut encore le débit thermique unitaire de l’eau chaude : 


19,8 


=p.q,.C..=10x x 4180 
dtc Pe dvc pe 3600 





4x=23.10 W/K 





d’où : 
NUT = 1120 x 230 
23.10 
NUT = 0,76 











5b). Efficacité 
Pour calculer E, il nous faut aussi le facteur de déséquilibre : 


R= dt min — ite 





dt max dy 
Les deux fluides ont la même masse volumique et la même chaleur massique, donc : 
Rp- 19,8 

dyf 82 
R = 0,24 











Les efficacités partielles E, et E, ont été déterminées dans la question 1 : 


1+2%X 0,24 2 


2 
Ë2 = 1=2%x 0.24 
1 — 0,48 exp - SS x 0.76) 
E, = 0,297 


On a montré que : 


E = 0,29 — 0,29 x 0,297 + 0,297 
E = 0,50 














Sc). Températures de sortie 

Toujours d’après la première question : 
Ty, = Te + RE(Tce — Te )= 20 + 0,24 x 0,50 x 70 
T = 28,4°C 














Tos = Ps Te) (Ei = E, E, + E; ) Siç E (Tee Tre) 
T,, E 90 — 0,50 x 70 
T..= 55°C 














On doit enfin revenir sur l’évaluation préalable des températures moyennes. Avec les 
résultats obtenus, ona: 


20 + 28,4 pans 
Tmf = Nr Që = 24,2°C pour une estimation de 25 °C : c’est très correct. 
0 + 55 ss 
Ts —— = 72,5°C pour une estimation de 70 °C : c’est assez près. 


Une itëration n’est donc pas indispensable. 


Commentaires 


Le problème combine un écoulement longitudinal dans un faisceau de tubes et un 
montage en réseau. 

Les deux demi-échangeurs ont des efficacités très voisines, malgré leur différence de 
catégorie (co-courant et contre-courant). Ceci est du au fait que le NUT est petit, ainsi que R 
(voir l’influence de ces deux paramètres dans les exemples des figures 3.1 et 3.2). On observe 
aussi que l’efficacité totale est nettement supérieure à chacune d’entre elles. 


PROBLÈME N° 14 : Réseau à 3 fluides 
Énoncé 


On considère un réseau ouvert mettant en œuvre deux échangeurs (fonctionnant à 
contre-courant) et trois fluides, avec un montage en série sur le fluide froid (fig. P14.1). 


Les températures d'entrée Ty, , Te, Te2 et les débits thermiques unitaires des trois 
fluides (qy, Qi, Qie2) Sont donnés, ainsi que les efficacités E; et E, des ëchangeurs. On 


impose également qy < Gjej €t Gif < qez- 


1. Calculer les températures Ty), Ty, Tes, Tes2 en fonction des données, puis le flux 


total échangé P. A partir de l'expression de Ð, préciser quelles conditions doivent satisfaire 
les températures d’entrée pour que l’ensemble fonctionne de façon correcte. 


Te 


BA. 





Fig. P13.1 


2. On se place dans le cas où q1 = q2. Les efficacitës des deux ëchangeurs sont 
notées E; et Ep. 

Comparer les puissances globales échangées Pap et Bpa selon que l'échangeur A 
est placé en position I et l'échangeur B en position 2, ou vice-versa (on étudiera le signe de 
la différence Pap — Pga). En déduire les règles d’une disposition optimale pour le réseau 


lorsque les températures d'entrëe sont imposées. 
3. Application pratique : 


Le fluide froid est de l'air, avec un débit q, = 800 m? /h. Les deux débits chauds 
sont égaux (cas de la question 2). On donne : Ty, = 25°C ; Th = 80°C ; Tha = 65°C. 

Les efficacités sont : E4 = 0,38 ; Ep = 0,25. 

Comparer les valeurs Pap et Ppa selon la position relative des deux ëchangeurs. 
Calculer les températures de sortie froides Ty( AB) et Ty( BA). Vérifier que les résultats 


sont en accord avec les conclusions de la question 2, et que la condition de fonctionnement 
(question 1) est satisfaite. 


Solution 


1. 
la). Sorties de l'échangeur N° 1 
Le débit thermique unitaire minimum dans l’échangeur est celui du fluide froid : 


dif < Gjej, d'où l'efficacité : 














: 1 pyt 
Dans ce cas, on a aussi d’après (3.11): Ep = T E., donc ici : 











Teer ~ Tos j ; 
Ej = R; E, e 0 ci avec R; = Armia = AP 
Ta E T to dt max dtcl 
If 
Ps F Teel E E, (Teel — T to ) 
Qici 











1b). Sorties de Vëchangeur N° 2 
Les conditions sont les mêmes : g,; < q,2. Le fluide froid entre à Tp; : 


E, = Epo = ——— 
2 Tez- Tri 

SI res Tr) 

En remplaçant Ty, par son expression obtenue en 7a) on obtient : 
TË ST pe FEST = Te) + EE; )( Ter Th) 














On a également d’après (3.11) : 














To = Ta i 
E. = R, E, = See Tes? avec R, = ay 
Less Ty; ic2 
dy 
Tez = Tat ES (Tes Tij) 
dic2 
En remplaçant Ty, : 
di 
Pit 7 E, Ps T T to Eten T to ) 
tc2 











Ic). Puissance 
Le plus simple est de calculer la puissance du système en faisant le bilan sur le fluide 
froid (circuit série) : 





D = qy [Ez (Toez -T pe ) + Ej (1 — Ez ) (Toer = Te] 








ld). Condition de fonctionnement 
Pour un fonctionnement correct, dans chaque échangeur, la température d’entrée 
chaude doit être supérieure à la température d’entrée froide. 


a TL 
Tez > Ty] 
Ceci est évident pour le premier. Pour le second, on devra donc vérifier que : 
Da > Tet ET. Fe) 


Les conditions de fonctionnement de l’échangeur N° 2 doivent être définies après 
celles de l’échangeur N° 1. 


2. 
Montage avec A en l et B en2: 


D’après l’expression de ® obtenue dans la question 1 : 
Pap “dy EN (Tao = T to JË E Te Ep ) (Toer E T to )| 


Montage avec Ben I et A en 2: 


Pra = dy (E (Toez = Te ) + Ep(1— E4 ) (Toer = The ) 
Après simplification, on obtient pour la différence : 











Pan PBA = qg (Ea Es (Teer Tea) 





E4 > Ep et Te > To 
Pap > Bgy Si et vice-versa. 
En < Ep et Taci < To 


On en dëduit une rëgle d'optimisation : 


- Si Toe] > Tea : placer l’échangeur le plus performant en position 1 


- Si Ty < Ta : placer l’échangeur le plus performant en position 2 
- Si Toe; = Leo Où E4 = Ep : Pap = Bpa. La disposition est indifférente. 


3 
II faut d’abord déterminer le débit thermique unitaire du fluide froid : 


If = Am Cf = PI af Cof 
La température moyenne de l’air est inconnue. Une valeur de la masse volumique 
Pf=110 kg/ m? (correspondant à T„ = 37°C ) est plausible pour démarrer le calcul. 


Pour la chaleur massique, peu dépendante de la température, on prendra 
Cf =1006 J/kg K. 


800 
ct x 1006 
Ki 3600 





dy = 246 W/K 


Disposition A-B 
D'aprës la question 2 : 
P ag = dif (Ep (Tee Es T to )+ Ea (1- Eg )(Teei — T to )| 
Pap = 246 [0,25 x (65 — 25 ) + 0,38 x( 1 — 0,25 ) ( 80 — 25 )] 





DAB =6315 W 











Du bilan sur le fluide froid : 
p = qy (Ts — Tre ) 
on tire : 





p 
Tp AB) =T p + —Ë zia ue 
qy 246 


Ty( AB) - 50,7 °C 














Disposition B-A 
Pa = qy LEA (Toez -Tje ) + Eg(1- En ) (Toei -Tp )] 
Pra = 246 x [0,38 x (65 — 25) + 0,25 x (1 -— 0,38 ) (80 — 25) | 














On calcule Ty, de la même façon : 


p 
TABA)E Te oe 000 
dy 246 


Ty( BA) = 48,7 °C 














On a ici TË.) > T2 : la meilleure disposition est bien celle où E, > E,, donc A-B. Le 
manque à gagner avec B-A dans ce cas particulier est de 485 W, soit 7,7 90. 


La condition de fonctionnement (question 1) est : T2 > Ty} 
Tti ST jet E) (Teej — Tye ) avec E) = E4 

Tsi = 25 +0,38 (80 — 25) 

T py = 45,9 °C 


La condition est vérifiée puisque 7.9 = 65 °C. 


Chaleur de vaporisation Ly de différents corps à la pression de 1 bar 


Désignation Ly (kJ/kg) 
Air 197 
Alcool éthylique 846 
Alcool méthylique 1101 
Ammoniac 1369 
Butane 402 

Eau 2256 
Ethane 490 
Ethanol 846 
Fluides frigorigènes non chlorés 120 à 240 
Fuel domestique (à 100 °C) 260 
Méthane 511 
Propane 448 
Toluène 356 


Tension superficielle de l’eau (à Ta = 373K ): © = 38.10 N/m 


Propriétés de la vapeur d’eau surchauffée 
Pression : 1 atmosphère (température de saturation = 373 K) 















0,216 x 10 
0, 242 
0,311 
0,386 
0,470 
0,566 
0,664 
0,772 
0,888 
1,020 


Lord 
400 | 0,554 |1,34 
450 | 0,490 |1,53 
500 | 0,441 |1,70 
550 | 0,400 |1,88 
600 (0,365 |2,07 
650 | 0,338 |2,25 
700 | 0,314 j2,43 
750 | 0,293 |2,60 
800 | 0,274 |2, 79 


































Propriétés thermophysiques de différents fluides (d’après J.Taine et J.P.Petit : 
Transferts Thermiques, Mécanique des fluides anisothermes, Dunod, 1989). 


Ke 4 bar 


a 






























































ae 1,13.107 1,57.1075 
1,85.10-5 Ú 57.1075 9 22105 
2,08.107$ | 2,08.1075 2,98.1075 
2, 29.1075 2 59.1075 3,76.1075 
2,48.1075 | 3,16.10-5 L,64.1075 
2,67.10-5 |3,79.10-5 5, 57.1075 
2,85.1075 | 4,43.1075 6,53.1075 

5,13.10-5 7,51.1075 





5,85.105 8,58.1075 
9,67.1075 
10,8.1075 
12,0.10-5 
13,1.1075 
14,3.1075 
15,5.1075 
16,8.1075 
19,5.1005 
22,0.1075 
24,8.1075 














0,757.1075 
1,563.10 75 


1,02.10-5 
2,21.10-5 
3, 74.1075 
5,53.1075 

















11,72.10-5 
13,60.10- 
15,61.107° 





0, 581.1075 
0, 832.107* 
1, 119.107 
1, 439.107" 
1, 790.1075 
2, 167.107 
2,574.107ë 
3, 002.1078 






4,48.1075 


T P p v 
K |kg/m?| kg/m.s m?/s 


HYDROGENE 


B 
Jika. K Ta K eag 1/Kelvin 
























































0,0981 | 7,92.10- | 8,06.10-5 0,156 
0,0819 | 8, 96.107€ | 10,9.10-5 0,182 
0,0702 )9,95.10-5 | 14,2.1075 0.206 
0,0614 | 10,9.10-5 | 17,7.10-5 0,229 
11,8.10-6 i 21,6.10-5 0,251 

25, 7.1055 0,272 

30,2.10-5 0,293 

35,0.10-5 0,315 

45,5.10-% 0,351 

56,9.10-5 0,384 

69,0.10-5 0,412 

82, 2.1075 0,440 


96, 5.1075 0,464 





























0,795.10-5 
1,144.1005 
1,586.1075 
2,080.10-5 
2,618.1075 
5 :3,199.10-5 
3,834.10-5 
4,505.1075 
6 )5,214.10-5 


0,0182 
0,0226 
0,0267 
0,0307 
0,0346 
0,0383 
0,0417 
0,0452 
0,0483 











ETHYLENE GIVCOL 
1130 | 65.1073 | 5,75.10-5 
1117 |21,4.10-% | 1,92.1075 
1101 j9,6.1073 10,869.10-5 
1088 | 5,2.107 )0,475.10-5 
1078 | 3,2.1073 10,298.10-5 
1059 12,15.10-% |0, 203.1075 













































899 3,850 | 4,28.107% 













888 | 0,80 0,90 

876 | 0,21 0,24 

864 | 0,072 | 0,0839 
852 | 0,032 | 0,0375 
840 | 0,0170 | 0,0202 
829 | 0,0102 | 0,0123 
817 | 0,0065 | 0,0080 
806 | 0,0045 | 0,0056 


siak (jat mije | 7 grain 


1002 i m 10-7? | 0,179.10- 4218 0,552 |113,1.107#1 13,6 |0,66.1074 
K 1001 1,30. 1073 | 0,130.105 4192 ‘| 0,586 |13,7.10-8 | 9,30 |0,88.10-* 
20 | 1001 | 1,00.10-5 0, 101.10-5 4182 0,597 14.3.1085 7,02 2. 06.1074 
40 | 994,6 )0,651.1031 0,0658.10”5 4178 | 0,628 |[15,1.10-8| 4,34 |3,72.107* 
60 | 985,4 |0,469.10-3| 0,0477.10-5 4184 0,651 |15,5.10 8 | 3,02 |5,15.107* 
80 | 974,1 |0,354.10-3|  0,0364.10-° 4196 0,668 |16,4.10-8| 2,22 |6,55.107* 
100 | 960,6 [0,281.10-#| 0,0294.10-5 4216 0,680 |16,8.10-% 1! 1,74 |7,49.107* 
120 | 945,3 10,234.10-51 0,0247.10-5 4250 0,685 |17,1.10-8 | 1,446 | 8, 92.1074 
140 | 928,3 |0,198.10-%| 0,0214.10-5 4283 0,684 |17,2.1078 | 1,241 | 10,0.107* 
160j 909,7 Jo,172.105) 0,0189.10-5 4342 0,680 117,3.10-% | 1,099 | 10, 7.104. 
180 | 889,0 j0,154.10”3j 0,0173.107 4417 0,675 |17,2.107ë | 1,004| 11, 4.1074 
200| 866,7 )0,138.10-5j 0,0160.10-5 4505 0,665 |17, 1.107ë | 0,937 | 14, 1.1074 
220| 842,4 )0,125.10”5) 0,0149.10-5 4610 0,653 116,8.10”5 | 0,891 15, 0.1074 
240| 815,7 |0,117.10-$| 0,0143.10-5 4756 0,635 116,4.10”8)0,871)18,0.107“ 
2601 785,9 j0,108.10”5) 0,0137.10-5 4949 0,611 |15,6.10-8 | 0,874 |21, 3.1074 
280| 752,5 |0,102.1073| 0,0135.10-5 5208 0,580 114,8.1078 | 0,910 | 26, 8.107“ 
300| 714,3 [0,096.10-%| 0,0135.10-5 5728 0,540 !13,2.1078 11,019 
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- à courants parallèles, 1.2 ; 4.3.2 ; problème 13. 

- à fluide isotherme, 2.4 ; 3.5.2 ; 4.2.1.2 ; 4.2.1.3 ; 
problèmes 4, 5, 10. 

- à modules, 1.3.1 ; 7.2.4.2. 


- à passes, 1.3.2 ; 3.7.14 ; 7.4.2 ; 
problèmes 8, 9, 10. 

- à plaques, 1.1.2 ; 1.3.2 ; 6.1 ; 7.4.3 ; problème 7. 

- à tubes ailetés, 4.3.4 ; problème 11. 

- bitubes, 4.4 ; 8.2 ; problème 3. 

- co-courant, 2.2 ; 3.4.1 ; problème 13. 

- compacts, 4.3.4.4 ; problème 11. 

- en épingle (à faisceau en U), 1.1 ; problème 13. 

- en série, 7.2 ; 7.4.1 ; problème 12. 

- en série-parallèle, 7.3 ; problème 12. 

- équivalents, 7.2.3 ; 7.2.5 ; 7.3.4 ; 8.6.1. 

- P-N, 1.3.2 ; 3.4.3 ; 3.7.2 ; 3.7.3 ; 7.2.4.1. 

- tubulaires, 1.1.1 ; 4.3.1 ; 4.3.3 ; 742; 
problèmes 6, 8, 9, 10, 13. 

Écoulement 
- autour de faisceaux de tubes, 4.14V ; 4.3.1. 
- de transition, 4.2.3. 


- laminaire, 4.2.1 ; 4.4.2 ; 5.2.2 ; 85. 
- turbulent, 4.2.2 ; 4.4.3 ; 8.4.2 ; 8.5. 


Efficacité 
- d’un échangeur, 3.2 ; 3.4.3 ; 7.2.3 ; 7.3.4. 
- d’une ailette, 4.3.4 
- moyenne en régime variable, 8.8. 
- relative, 3.2 ; 7.2.2 ; 8.8. 
Encrassement, 6.1 ; 6.2.3 ; 8.8. 
Évaporateurs, 2.4 ; 5.3. 


F 


Facteur de déséquilibre, 3.2 ; 7.2.1 ; 8.8.2. 

Facteur de forme, 4.4.1. 

Fluide brassë, 1.2 ; 3.4.3 ; 3.5.3 ; problème 8. 
Formule de Colebrook, 4.2.2.3 ; problème 9. 


L 


Longueur caractéristique, 4.3.1.5 ; 4.3.2 ; 4.5. 
Longueur de référence, 4.1.5 ; 4.2.2.6 ; 4.3.1.5 ; 4.5. 


M 
Méthode NUT, 3.6 ; 3.7 ; 7.2 ; problèmes. 
N 


Nombres de Nusselt, Péclet, Prandtl, Reynolds, 
Stanton, 4.1.5. 
NUT, 3.3 ; 3.4.3 ; 6.2 ; 7.2.5 ; 7.34. 


P 


Pas d’un faisceau, 4.3.1.1 ; problèmes 9, 13. 
Perte de charge, 4.2.2.3 ; 4.3.2 ; problèmes 9, 10. 
Pincement, 3.7.1 ; 3.7.3 ; 7.6. 

Puissance d’un échangeur, 3.1 ; 6.2. 


R 


Réseau, problèmes 12, 13, 14. 
- à fluide intermédiaire, 7.5.1. 
- maillé, 7.4. 
- ouvert, 7.5.2. 

Rugosité, 4.2.2.3. 


Température 
- de film, 4.1.3 ; 4.1.4 ; 5.1 ; problème 10. 
- de mélange, 2.1 ; 2.2.9 ; 4.1.2 ; 4.14. 
- de référence, 4.1 ; problèmes 8, 9. 
Temps de retard, 8.3.1 ; 8.6.3. 


Thermodépendance, 4.2.1.1 ; 4.2.2.1. V 
Tubes 


- corruguës, 4.2.2.5 ; 4.4.4. Viscositë, 4.2.1.1. 
- faisceaux, 4.3.1.1; 4.3.1.6; 4.3.3; 43.43; Vitesse massique, 4.1.5 ; 4.3.1.2. 
problèmes 9, 13. Vitesse de référence, 4.3.1.2 ; 4.3.4.3. 


- nappes, 4.3.1.1 ; problème 9. 

- section annulaire, 4.4 ; problème 3. 

- section circulaire, 4.2.1.2 ; 4.2.2.2. 

- section rectangulaire, 4.2.1.3 ; problème 9. 


